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ABSTRACT 


Radio wave propagation in the earth's crust is reviewed. This has been 
studied by many investigators who usually assume the crust to be a planar 
slab of low loss-tangent, bounded on top and bottom by regions of large loss- 
tangent. Available conceptual depth profiles of complex dielectric constant 
of the lithospheric waveguide are reviewed and theories extant summarized. 
Sources were usually VED's inserted through vertical bore-holes. For ocean 
emplaced antennas located near the ocean-lithosphere boundary, the sources 
are horizontally polarized and propagation is predominantly by means of lateral 
waves. Calculations are reported for a wide range of conditions. The useful 
characteristics of lateral wave transmission are discussed. The measured 
properties of scaled waveguide propagation that have been reported could not 
accurately scale the losses in the overburden since it was represented by a 
metal plate. The principles of electrodynamic similitude are delineated and 
two different experiments using similitude with different depth profiles are 
discussed, t 

To study the effects of boundaries on transmission (including wall losses) 
and simultaneously to make measurements of embedding and moving antennas in 
them, the laboratory models should have "walls" composed of a liquid with 
suitable excitation. For an orderly build-up of modelling measurements, a 
laboratory set-up is described. Included would be measurements with boundaries 
that are sloping, undulating or irregular. Though available materials may 
limit accurate modelling according to similitude scaling, meaningful and 
illuminating information could be obtained for better understanding of 
lithospheric transmission. The more complex depth profiles require a much 
more elaborate set-up, patterned possibly after that of lizuka. 
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zontally polarized and propagation is predominantly by means of lateral waves. 
Calculations are reported for a wide range of conditions. The useful charac¬ 
teristics of lateral wave transmission are discussed. The measured properties 
of scaled waveguide propagation that have been reported.could not accurately 
scale the losses in the overburden since it was represented by a metal plate. 
The principles of electrodynamic similitude are delineated and two different 
experiments using similitude with different depth profiles are discussed. 

To study the effects of boundaries on transmission (including wall losses) and 
simultaneously to make measurements of embedding and moving antennas in them, 
the laboratory models should have "walls" composed of a liquid with suitable 
excitation. For an orderly build-up of modelling measurements, a laboratory 
set-up is described. Included would be measurements with boundaries that are 
sloping, undulating or irregular. Though available materials may limit accu¬ 
rate modelling according to similitude scaling, meaningful and illuminating 
information could be obtained for better understanding of lithospheric 
transmission. The more complex depth profiles require a much more elaborate 
set-up, patterned possibly after that of Iizuka. 
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1. Introduction 

Interest in the electrical and geophysical properties of earth's crust 
has been recently encouraged by the possible existence of a lithospheric 
plane waveguide with conductivities lower than previously accepted. The 
region of the earth's crust forming this low conductivity waveguide was 
assumed to be a planar slab, 10 to 30 km thick, bounded on top by an over¬ 
burden (or ocean) of high conductivity and on the bottom ("underburden") by 
a high conductivity mantle or Moho (as it is often referred to in the 
literature). A review of the crustal properties was given at an ONR 
Lithospheric Workshop in mid-March, 1978, chaired by J. Heacock with papers 
presented by G. Keller, G. Simmons, P. R. Bannister, and C. Cox amongst 
others. 

There have been some excellent reviews of the field in recent years. 

For example. Wait (1971) in the AGU Monograph No. 14, covered principally 
the waveguide analytical approaches. And 3annister (1978) gave an excellent 
review of the idealized waveguide (homogeneous profile) of Wait and Spies 
plus the results of Field and Dore for non-homogeneous profiles of Levin, 
Keller, and Housley. We wish merely to supplement here some of their 
writings. 

While most of the analytical work tended to treat waveguide modal 
solutions, there exists the possibility of lateral wave transmission between 
ocean emplaced HED transmitting and receiving antennas located near the 
ocean-lithosphere boundary. 

It was desired to study the feasibility, in laboratory models, of 


measuring the properties of lithospheric transmission, including wall losses, 
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using the principles of electrodynamic similitude for scaling. Scaled 
models have been used for a homogeneous waveguide in one case, and in another 
for a non-homogeneous profile of conductivity increasing with depth. Both 
cases used metallic upper plate upper walls which do not accurately scale 
the properties of the overburden or ocean. There is a question of available 
materials for such modelling. 

In Chapter 2, we present separately the various available conceptual 
depth profiles, starting from the 1961 work of Wheeler and that of Mott and 
Biggs in 1963, both works being expositions of concepts of communication 
from land via "A deep waveguide" in the crust to deeply submerged submarines. 
There follows a series of descriptions of several non-homogenecus depth 
profiles. 

Chapter 3 contains a summary review of the various theoretical works on 
the profiles of Chapter 2. Most of these concern modal solutions for wave¬ 
guides with plane boundaries; results are often illustrated with boundaries 
having infinite conductivities. In some cases, the lower boundary is sharp; 
in others the complex dielectric constant is tapered to increase with in¬ 
creasing depth. Finally, the work with antennas emplaced in the ocean near 
the ocean-lithosphere boundary is cited by way of introducing the next 
section. 

In Chapter 4 is an expose with numerical calculations of lateral waves 
from an HED source emplaced in the water near the ocean-lithosphere boundary. 

A previous development of King and Sandler for an HED in water below air, 
which extended the work of Banos, is extended to the case where the lithosphere 
replaces the air. A wide range of parameters were used for calculating the 


fields 
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Chapter 5 contains a brief development of expressions for electrodynamic 
similitude. Examples are given for laboratory frequency and distance scaling 
by Brown and Gangi for a homogeneous depth profile of dielectric constant and 
conductivity, plus some results of lizuka with the non-homogeneous profiles 
achieved with an agar-agar gel diffused with sodium chloride. 

In Chapter 6 is contained a discussion of lateral wave modelling in 
combination with measurement of losses in the walls for this and for wave¬ 
guide excitation. A suggested laboratory set-up could be used for sloping 
boundaries or boundaries that are rough or undulating. 

A sunmary and conclusion with recommendations are contained in 
Chapter 7. 

There follows the list of references and supplemental bibliography in 


Chapter 8 
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2. Lithospheric Profiles 


2.1 Introduction 

Most of the earlier lithospheric waveguide depth profiles were based 
on the laboratory measurements of samples of rock as functions of tempera¬ 
ture and pressure. With few exceptions, emphasis was placed on the depth 
variation of conductivity o; this was due probably to the belief that the 
expected exponential attenuation would be low enough only at low frequencies 
where the loss tangent p = 60 c X Q /e r is quite large (a in Si/m, \ Q the free 
space wavelength in meters, and e r the relative dielectric constant). In 
analyses of these profiles, most investigators had assumed e r to be constant 
with depth in the lithospheric waveguide, with values of 4 to 10. 

Earlier analyses assumed the waveguide to have a depth of 10 to 40 km. 
On top, it was bounded by an overburden or sea of such high loss-tangent 
that this region was represented by a highly conducting metallic plate in 
many analyses and models. 

At the bottom, the earlier analyses assumed the conductivity towards 
the hot "Moho" discontinuity to be so high that it again was represented by 
a highly conducting plate. These earlier analyses assumed a homogeneous 
parallel plate slab, somewhat idealized, to calculate the fields in the 
guide. Refinements grew as more complicated profiles of the crust developed. 
We trace here the profiles from the early one of Wheeler to that of Wait, 
then to Levin and coworkers, to Keller, to Housley and to the work of 


Simmons. 
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2.2 Wheeler’s Profile 

Wheeler's (1961) early paper, interestingly enough, concerned itself 
with a means of communicating to a submarine near an ocean floor from a 
shore based source, A sketch of the concept is shown in Figure 2.1 (a); 
the transmitter, shown in Figure 2.1 (b), passes through about 1 km of 
overburden and then extends 2 km in a borehole into the waveguide where 
it excites a TEM mode. The conductivity profile is shown in Figure 2.1 (c) 
as a function of depth and temperature. The boundaries are roughly 2 and 
20 km deep, where the rock conductivity may be 10“ 6 to lO”^ 1 Si/m, a useful 
value being 10 -8 Si/m. The relative dielectric constant was 6. The "upper 
wall", perhaps 1 to several kms thick has a range of 4 Si/m for sea water 
down to 10 -4 Si/m in rather dry minerals. The lower region was termed the 
"inverted ionosphere" or "thermal ionosphere" by Wheeler (1961) with tem¬ 
peratures of 300° to 600° C. 

2.3 Mott and Biggs' Profile 

A couple of years later, Mott and Biggs (1963) also considered VLF 
propagation below the sea. Their concept is shown in Figure 2.2 (a) for 
the coiranunicating region, the upper boundary 1 km thick being the wet, 
highly conducting continental surface or the ocean. The propagating region 
was assumed to consist of dry, low-temperature low conductivity granites or 
basaltic rocks, 20 to 70 km thick. The bottom was formed by ultra-basic 
rooks, having higher conductivity because of the high temperature. 

Propagation takes place between 2 vertical electric dipoles (VED) 
extending 1 or 2 km into the low conductivity region, with an upper value 
of 10“ 6 Si/m for long-distance communication. 
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Earth conductivity estimates are shown in Figure 2.2 (b). The curve 
labeled number "1" indicates measured earth conductivities of 10 -3 for 
poorly conducting earth to 4 Si/m for sea water. Curve "2" is from Hughes 
(1953) for olivines, with McDonald's (1959) estimates of temperature 
distribution. Curve "3" is an estimate of maximum conductivity, by Garland 
and Webster (1960) for western Canada pre-Cambrian basement rocks. Curve "4" 
is that of Lahiri and Price (1939) while curve number "5" is due to Cantwell 
and Madden (1960). Some of these curves are included in a composite one 
drawn by Iizuka (1968), discussed later on. The final conductivity-depth 
profile assumed by Mott and Biggs (1963) is shown in Figure 2.2 (c). In 
essence this is a two-region model, the bottom being 10 - 20 km deep with 
a bottom region having a conductivity 10”® to 10 -7 Si/m, and a dielectric 
constant of 6 with a top region of sea water with 4 Si/m for conductivity, 
perhaps 1 to 2 km thick. 

2.4 The Profile of Brown and Gangi 

Brown and Gangi (1962, 1963) and Gangi (1966) gave a diagram of earth 
conductivity as a function of frequency and depth, as shewn in Figure 2.3. 

The surface and near surface overburden had conductivities ranging 
from 10" 1 to 10” 4 Si/m, with an average of 10” 2 to 10” 3 Si/m. The relative 
dielectric constant had a wide range of values, with a typical value of 4 to 
10 being widely used. 

For the crust, the values were inferred from tables in Von Hippel 
(1954) for the effects of temperature and water content of quartz and sili¬ 


cates (pages 303-314 of his book). The low conductivity value of granite 
is frequency sensitive and the authors show values for 1, 10 and 100 kHz. 














At 30 km depth, the conductivity increases rapidly, temperature being about 
300° C, relatively independent of frequency. Near the Moho the conductivity 
is about 4 x 10 -3 Si/m. Relative dielectric constant was of the order of 5, 
down to 30 km. 

In their theoretical studies, the "real earth" situation was approximated 
by a uniform slab between two highly conducting planes (See next Chapter). 

The conductivities chosen were according to frequency in the low conductivity 
region. 

2.5 Levin's Profile 

An idealized cross-section of the earth's crust was delineated and shown 
by Levin (1966) in Figure 2.4 illustrating the crust under a continent vs. an 
ocean basin, as well as the continental shelf. Under the continent is granite 
underlain by basalt down to the Moho. Proceeding through the shelf, the 
granitic rock disappears leaving basalt under the ocean and volcanoes. The 
Moho may be 35 km deep under the continent and about 10 km under the ocean. 

After a considerable discussion.. Levin (1966) showed a trial profile of 
both conductivity and relative dielectric constant with depth. Levin (1966) 
cited the measurements of rock characteristics by Acker and Mueller (1966) 
from 3 km drill holes associated with a full scale propagation experiment. 

His conductivity profile is summarized in Levin (1971) and shown in 
Figure 2.5. The variation of dielectric constant is discussed later in 
the paper by Schwering, Peterson, and Levin (1968), and is overlooked by 
many authors. 

In his AGARD paper, Levin (1966) concluded that "full scale studies of 
lithospheric propagation must be carried out in boreholes nearly 10 km deep 
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to assure access to the dielectric zone (loss-tangent less than unity) of 
the crust over the entire length of the path". 

In Figure 2.5, two interpretations are shown for conductivity vs. 
depth. The curve labeled "dry" is more nearly in favor in view of recent 
consideration of laboratory measurements and interpretations, especially 
those regarding the influence of microcracks (Simmons, 1978). 

2.6 The Schwering, Peterson, and Levin Model - Inverse Square Profile 

Schwering, Peterson, and Levin (1968) analyzed lithospheric propaga¬ 
tion based largely on the Levin (1966) model, the dry portion of conductivity 
of Levin (1971) shown in Figure 2.5. The variation of dielectric constant 
and conductivity are shown in Figure 2.6 as solid lines as the "experimental 
values"; the dotted curves are used for their analysis, described in the next 
Chapter. 

It was convenient to distinguish three zones in the crust, with elec¬ 
trical constants given for 1 kHz, and shown in Table 2.1, as well as other 
postulated characteristics. The waveguide is Zone 2 in essence. 

2.7 Wait's Exponential Profile 

After studying the performance of electrical constants in a number of 
ways. Wait (1966c), offered a unique profile suitable for analysis (Wait 
1966a, 1971), to attempt a quantitative result. His model is illustrated 
in Figure 2.7. In a cylindrical coordinate system, the plane z = 0 sepa¬ 
rates the highly conducting overburden having a complex relative dielectric 
constant K g , the overburden top being at z = -h. The basement has a rela¬ 
tive dielectric constant K (z) which is complex, having a constant relative 
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Table 2.1 Descriptions of crustal zones according to Schwering, Petersen, 
and Levin given in their 1968 article, Schwering, Petersen, and 
Levin (1968) 
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dielectric constant and a conductivity which increases exponentially with 
depth. The reference depth was chosen for convenience as that depth z d 
where the loss tangent was two; for VZF it appears that z^ is of the order 
of 30 km. At very great depths, where the loss tangent greatly exceeds 
unity, Y.(z) varies exponentially with depth. Near z * 0, the ideal case 
would have K(0) real. 

The results of using this profile will be given in the next Chapter 
where he compares his results with those of other authors, 

2.3 Iizuka's Profiles 

Iizuka performed a number of experiments using a "typical" lithospheric 
profile deduced from a number of sources. The typical case of a profile of 
loss-tangent with constant relative dielectric constant is shown in Chapter 
5, Figure 5.6. for a dielectric constant of 4 and a typical frequency of 
5 kHz. The compilation of conductivity profiles of Iizuka (1968 ) is shown 
in Figure 2.8, They were deduced from Wait (1963), Wheeler (1961), Mott 
and Biggs (1963), Keller (1963), Ames, Frazier, and Orange (1963), Cantwell 
and Madden (1960), Watt, Mathews, and Maxwell (1963), deBettencourt and 
Frazier (1963), Garland and Webster (1960) among others. 

In Figure 2.8, the conductivity is plotted as a function of depth to 
depths of 40 km or slightly more. The conductivity ranges from 10 -6 to 10 -1 
Si/m, the larger values being near the earth's surface or at great depths. 
Since he preferred an analysis based upon loss-tangent, he noted "lower 
limit of tan 5" near a conductivity of 10“ 8 Si/m and a higher limit near 
10 km depth with a conductivity of about 10“ 4 Si/m. (For his "typical" 
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loss tangent curve, Iizuka (1968 5 assumed a conductivity of slightly less 
than 10 7 Si/m with a dielectric constant of 4, at a frequency of 5 kHz) 

2.9 Keller's Profile 

Keller (1963, 1971, 1978a,b) has studied in depth the electrical crustal 
properties. His most recent profile of conductivity is shown in Figure 2.9. 

It has been widely used and quoted. It appeared first in Gallawa and Haidle 
(1972), as a private communication, again in Field and Dore (1973), discussed 
by Keller (1978a) at the ONR Lithospheric Workshop and is reprinted here from 
Keller (1978b). 

While the value of 10“ 4 Si/m at very shallow depths appears low (Field 
and Dore, 1973; Bannister, 1978), the curve down to about 10 km is based on 
Appalachian well-log data and an extrapolation thereof towards minimum values 
of the order of 10” 7 Si/m at 10 km. The lower portion has increasing values 
from those of granite to underlying basalt. A rather complete discussion is 
given in Field and Dore (1973). 

If water saturated microcracks exist in rocks at 5 to 8 km depth, mini¬ 
mum conductivities must be greater than 10 -5 to 10 -6 Si/m, and practical long 
range communication does not appear feasible. However, recent data concerning 
microcracks by Simmons (1978) tend to indicate conductivities could be lower, 
say 10 8 to 10 Si/m, than that previously believed (10 _/ Si/m). 

2.10 Housley's Profiles 

Housley (1973) made laboratory measurements on the conductivity of dry 
rock, see Simmons (1978). Results show even lower conductivities than those 
of Figure 2.5 of Levin and Figure 2.9 of Keller. Using the new conductivity 
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Figure 2.9 Keller’s profile of conductivity and temperature vs. depth. 

From Field and Dore (1973), first as a private communication 
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See also Bannister (1978). 
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data, the profiles are similar to those of Keller in Figure 2.9, for depth 

less than 10 km, and are shown in Figure 2.10. At greater depths, much 

_8 

lower conductivities are indicated (less than 10 Si/m). 

The three profiles correspond to geotherms in three different heat 
flow provinces (i.e., BR = Basin and Range; EUS = Eastern United States; 

SN = Sierra Nevada). The thermal gradients for the Keller model was 16°/km, 
while for the Housley provinces they were BR (24°/km), EUS (14°/km), and 
SN (9°/km), where the temperature is in (°C). 

2.11 Comment 

There are many differing profiles of conductivity vs. depth which are 
not necessarily definitive. Continued experimental work remains to be done 
in the laboratory - see Simmons (1978) - and in the field - see Keller (1978a). 



3. Summary Review of Theories 


3.1 Introduction 

When viewing the number of more recent depth profiles of electrical 
constants in Chapter 2, such as those of Keller in Figure 2.9 and of Housley 
in Figure 2.10, interest has been revived in the possibility of long-distance 
communication via a sub-surface waveguide. The minimum conductivities are 
several orders of magnitude lower than those postulated in the early 1960 
era. The lower 'onductivity values have received support from studies of 
the effects of microcracks by Simmons (1978). 

A number of review papers (Rjanzcev and Shabel'nikov, 1966; deBetten- 
court, 1966; Levin, 1966, 1971; Wait, 1966c, 1971; Gabillard, 1966; Gabillard, 
Degauque, and Wait, 1971; and Bannister, 1978) have been published on the 
various mechanisms dealing with lithospheric propagation, including ground 
waves, "up-over-and-down," and sub-surface waveguide mechanisms plus com¬ 
munication parameters of antennas, noise, etc. Here we shall briefly sum¬ 
marize those papers dealing principally with the crustal waveguide mechanism. 
Many focus attention on attenuation of the modes, a few give examples of 
field strength variations. 

3.2 Some Early Works 

It is interesting to realize that two of the earliest papers considered 
the possibility of using the lower-loss region of the crust for the purpose 
of illustrating a potential mechanism for communication with submarines at 


great depth 






Wheeler's (1961) very early concept was shown in Figure 2.1. He 
imagined the existence of a very deep waveguide. With a relative dielec- 
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trie constant of 6 and a conductivity of 10 to 10 Si/m, he obtained 
a total attenuation of the TEM launched wave of 0.067 dB/km. 


Mott and Biggs (1963) also considered undersea communication between 

two vertical electrical dipoles (VED) both inserted into lower conductivity 

crustal rock. They envisioned propagation at frequencies of the order of 

—6 —7 

100 kHz to take place in a plane half space of conductivity 10 to 10 Si/m 
bounded on top by the layer (ocean) having a conductivity of 4 Si/m. They go 
on to assume the lower region attenuates the wave but does not affect the 
propagation mechanism - surface wave along the interface like normal surface 
wave propagation in air above the sea but where the air has been replaced by 
the rock. See also Carolan and deBettencourt (1963). At lower frequencies, 
less than 10 kHz for example, the lower region behaves as a waveguide, and 
"that a TEM waveguide mode will be important", as Wheeler (1961) had 
postulated. Their conductivity profile is shown in Figure 2.2(c). 


3.3 Slabs with Uniform Depth Profiles 


For analytical studies of their depth profile shewn m Figure 2.3, 
Brown and Gangi (1962, 1963) and Gangi (1966) assumed the crustal waveguide 
to be homogeneous slab bounded by two half-spaces which were highly 
conducting.* For their analysis, the boundaries were assumed to have 
effectively infinite conductivity. The waveguide was excited by a vertical 


electric dipole.A mode series representation was used to obtain the field. 


♦Reference may be made to Figure 5.1 for a sketch of the guide 
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Calculations were made of the vertical electric field strength versus 
range for several low frequencies for a guide depth H = 35 km, a relative 
dielectric constant of 4 and zero antenna depths. A typical result is 
shown in Figure 3.1 for 1, 3, and 10 kHz, the conductivity values corres¬ 
ponding to their profile shown in Figure 2.3. For this configuration, the 
modal cut-off frequency 


f (m) = m v 2 /2H (3.1) 

where v^ is the phase velocity in the crustal waveguide and where m is the 
mode number, whence f (m) = 2.1m, kHz. Hence, for 1 kHz, only the TEM 
mode (m = 0) propagates readily; for 3 kHz there are two modes and for 10 
kHz, there are five modes (m = 0, 1, 2, 3, and 4), propagating with low 
attenuation. Interferences occur because each mode has a different attenua¬ 
tion and phase shift. 

Using Wait's (1957) mode theory for the earth-ionosphere waveguide 
with lossy ground, but replacing the air by a low-loss waveguide, deBetten- 
court (1962, 1966) calculated the vertical electric field versus range for 
infinitely conducting walls for a guide width of 20 km, relative dielectric 
constant of 4, and for conductivities of 0 and 10 6 Si/m; the dipole moment 
was 16,900 amp.m. (A plot of the results may be seen in Figure 6 of Wait, 
1966c) . 

Spies and Wait (1971) and Wait and Spies (1971) have extensively con¬ 
sidered the wave in an idealized waveguide. From a general solution of 
multiple horizontal stratified layers (Wait, 1962), Spies and Wait (1971) 
chose a restricted 3 region case, shown in Figure 3.2, in which the 
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"overburden" region had infinite conductivity, the crustal "waveguide" 
(region 2 ) had uniform conductivity of c 2 = a and relative dielectric 
constant K 2 , and the "sub-basement" (region 3) had a conductivity O 3 = a c 
and relative dielectric constant of K c , and was infinitively thick. The 
waveguide has a thickness h 2 = h = H. The boundaries of the various 
regions are sharp. In their calculations using a Newton method of suc¬ 
cessive approximations, Spies and Wait (1971) and Wait and Spies (1971) 
considered frequencies from 1 to 100 kHz; these were expended later for 
frequencies down to 1 Hz (Wait and Spies, 1972). 

Wait (1966, 1971) demonstrated the principle of "super-position of 
losses," which also was used by Watt, Leydorf, and Smith (1966). For 
modes of sufficiently low attenuation, the total exponential attenuation 
of a mode can be considered as the sum of three terms, i.e», for a mode m 

A t (m) = Ap + + Aj (3.2) 

where the three components are respectively the plane wave, upper boundary, 
and lower boundary contributions to the total attenuation rate. 

We are generally concerned with a crustal waveguide having a low- 
loss-tangent region whence it can be shown that 

A p = 1.64 x 10 6 n 2 /JT 2 dB/Kn (3.3) 

and is independent of frequency if the electrical "constants" are constant. 
For a given plane slab boundary, the wall attenuation contribution depends 
upon the guide loss-tangent and whether the boundary is electrically thick 
or thin. Expressions based on approximation similar to those of Wait (1963) 
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have been given for several cases by Watt et al. (1966) . As aforementioned, 
we are generally interested in the crustal waveguide which has a low loss- 
tangent, and the boundary layers have a high loss-tangent, and are electri¬ 
cally thick. For this case, if o is the conductivity of one of the uniform 
slab boundary layers, then Watt et al (1966) give the boundary attenuation 
for the TEM (m = 0) mode as 

A 0 (0) = 2.29 x 10“ 2 (K 2 f/a )Vh dB/Km (3.4) 

Whilst for the TMj Mode (m = 1) 

-‘h 

A b (1) = 4.58 x 10" 2 (X 2 f/° Q ) h [l-( f c / f)2 J /H dB/Km ( 3 ' 5 > 

where f is f (m = 1) from equation (3.1) 
c c 

In Figures 3.3 and 3.4, taken from Spies and Wait (1971), the total 
attenuation for modes of lowest attenuation are shown as functions of 
frequency for a guide thickness of H = 20 km, and upper layer of infinite 
conductivity, the waveguide with a dielectric constant K = 9 and conduc¬ 
tivities of 10 -6 or 10~ 2 Si/m, and the sub-basement layer with a uniform 

dielectric constant of K =15 and a conductivity of 10“ 2 Si/m (Fiaure 3.3) 

c 

or 10“ 3 Si/m (Figure 3.4). 

In a manner similar to that of Bannister (1978), we made calculations 
for a frequency of 10 kHz, the model cut-off frequency being 2.5m kHz. 

(The loss-tangent of the waveguide region is very small. The plane wave 
attenuation A^ for = 10“^ Si/m, = 4 works out to be .0546 dB/Km, 
whilst for 0 = 10" 6 , it becomes 0.546 dB/Km. The plane wave attenuation 

is critically dependent on conductivity. 













Frequency (kHz) 


Fig.3.4 Sane as Fig.3.3, but for (T. » 10 -3 Si/m (^pies and '’ait,l°71) 
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For the Spies and Wait (1971) model calculations for Figure 3.2, the 

upper layer conductivity is infinite so A (m) = 0 and A (m) = A ♦ A s (m). 

u c pi 

For H =* 20 km, Kj = 9, for the mode (m = 1) , f » 2.5 kHz. For 
= 10 ~ 2 Si/m and f = 10 kHz, the boundary attenuation is .007 dB/Kta 
while for a g ~ ^2 = Si/m it is .023 dB/Km. The higher the conductivity, 

the lower the attenuation or leakage through the wall. 

Similar calculations were made for the TEM (m = 0) mode. Some of the 
results are shown tabulated in Table 3.1 where the values in parentheses 
are rounded off tabular values taken from Spies and Wait (1971); the approxi¬ 
mate attenuation A T (0) and A^, (1) can be compared with Figures 3.3 and 3.4; 
the agreement is quite good here. 

It is interesting to consider one aspect of the attenuation of the 
lowest order mode, the TEX mode (m = 0), with reference to equations (3.2), 
(3.3), and (3.4). If the homogeneous slab of relative dielectric constant 
K., and conductivity C., has a thickness of 2h_ = 2h * 2H, but is bounded on 
top and bottom by homogeneous half-spaces of the same conductivity o a and 
dielectric constant the total attenuation (0) will be unchanged 

from that derived for the configuration of Figure 3.1. This will be referred 
to in Chapter 6 in the discussion of model measurements of wall losses with 
an image plane for region 1. This point has also been noted by Wait and 
Spies (1972) in their discussion of ELF crustal waveguide attenuations for 
the lowest order mode. 

The model just discussed is mindful of the analysis for a three layer 
structure, shown in Figure 3.5(a), as considered by Viggh (1963). His slab 
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TAPLF 3.1 


Aoproximate versus exact modal attenuations for modes of lew attenuation, for 

for idealized wavecruide shewn in Fiaure 3.2 
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10- 6 

.546 

.00^ .55 (.55) 

10-6 

.546 

.007 

.55 (.57) 

10“ 7 

.055 

.004 .059 (.059) 

10- 7 

.055 

.007 

.062 (.064) 


°~3 = 10-3 cj/p 

10-6 .546 .011 .56 (.56) 10~ f .546 .023 .57 (.59) 

10- 7 .055 .011 .066 (.071) l"" 7 .055 ."23 .079 (.090) 


For coroarison, the ficures in oarentheses are from Fnies and T- T ait,l°71, tables. 
The attenuations A are in dP/km. 












3.12 


of depth D with conductivity c 2 = 10“ 6 Si/m was bounded by identical top and 
bottom half spaces with conductivities = 10“ 2 Si/m; all three regions 

had the same relative dielectric constant of 10. The slab was electrically 
thin so that only the TEM mode was considered dominant. An example of the 
attenuation as a function of thickness is shown in Figure 3.5(b) where fre¬ 
quency is the parameter; both scales are logarithmic. 

Waveguides having low loss-tangents but bounded by layers on top and 
bottom have been considered by Gabillard (1966) and in a review paper given 
by Gabillard, Degauque, and Wait (1972). Essentially these papers refer to 
the work in France where they were able to demonstrate the propagation 
effects in mine galleries and in a uniform gypsum layer. 


3.4 

3.4.1 

Besides the aforementioned works of Wheeler (1961) and Mott and Biggs 
(1963) another early paper was that of Wait (1963, 1971)who considered the 
analysis of the exponential depth profile shown in Figure 2.7 in order to 
obtain quantitative results. For z positive, the complex relative dielec¬ 
tric constant was chosen for convenience to be exponential in form and 
given by 

K(z) = K q £l-}2 exp B(z-z d )j (3.6) 

where z^ is the waveguide "depth" which is discussed below. 

For z = 0, at the top of the waveguide 

K(0) * K Q [l-j2 exp (-Bz d )] * Kq 



(3.7) 


(a) Three-layer structure for analvsis 



Fig.3.5 (a) Three-laver structure for analvsis: (b) Attenuation =■< for electrical 

thin waveguide layer vs. laver thickness D , with freauencv as the 
parameter. (Viggh,1963) 











provided Bz d exceeds 3. With z d of the order of 30 km and B of the order 
of 0 . 2 /km, this restriction is true, i.e., the loss-tangent is extremely 
small and K 0 is essentially a real quantity at z = 0 . 

For z = z d , 

K(z) = K 0 (l-j2) = K q /5 exp (3.8) 

That is, the loss tangent at z = z d is 2. 

For very large depths 

K(z) = -jjK(z)j = -j2K Q exp B(z-z d ) (3.9) 

Since K q is real, the conductivity may be expressed as a function of 
depth as 

a(z) = K q exp (z-z^J/30 X Q (Si/m) (3.10) 

where X 0 is the wavelength in free-space in meters. These expressions 
assume that this real relative dielectric constant and conductivity are 
not functions of frequency. 

At z = 0, where the loss-tangent of the guide is very small, we desig 
nate the wavelength as 

\' = X /K 1 * (3.11) 

2 o o 

Wait (1963) derives the vertical field strength for a VED in the wave 
guide from a vertically directed Hertz vector. From the residues at the 
poles of the denominator in the expression for the latter, he obtains the 

i 

♦Wait uses the quantity X 0 to designate what we have called X 2 
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model resonance condition, finally expressing the field E, as E w where E 
is a reference field. This reference field is that of a vertical dipole 
in a dielectric half-space erected on a perfectly conducting ground plane 
of infinite extent. The expression for w is dimensionless and contains all 
the characteristics of the waveguide modes. Because the lower boundary is 
not sharply defined, the mode of lowest attenuation is not the TEM (m = 0) 
mode but one of the TM modes (m = 1, 2 . . .). 

For purposes of illustration. Wait (1963, 1971) calculated the attenua¬ 
tions for modes m = 1 and 2 for the special case where the overburden was 

so largely conducting that was assumed to be infinite. He chose the 

/ 

case where X ^ ~ 15 km. The attenuations A^ and A^ for the nodes of lowest 
attenuation are shown plotted as functions of the quantity B in Figure 3.6(a) 
and (b), respectively. The parameter on the curves is the waveguide depth 

i 

in "wavelengths" denoted by z^/X.. 

Wait (1963, 1971) shows how these curves can be scaled for other wave¬ 
lengths, for a given guide thickness in wavelengths. In these articles and 
another for plane slabs (Wait, 1966), he illustrates how the total attenua¬ 
tion for the lowest order modes may be shown to be the sum of the three terms, 
the plane wave attenuation plus the losses in the two "walls" as given in 
equation (3.2), (Wait, 1971). 

For the dominant mode attenuation, he notes the comparison of values of 
Schwering, Peterson, and Levin (1968) for their inverse square profile of 
0.256 dB/km compared with that value (0.275 dB/km) they obtained using Wait's 
(1963) method. Wait's (1963) values, he notes, lead to lower attenuations 


than those of Schwering et al (1968) because of the steeper gradients. 
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3.4.2 The Inverse Square Profile of Schwering, Peterson, and Levin 

The depth profile of Schwering, Peterson, and Levin (1968) is shown in 
Figure 2.6, the dashed curves being used for analysis. Both the real rela¬ 
tive dielectric constant K and conductivity a vary as the inverse square of 
the depth z measured positively downward. The vertical electric dipole 
(VED) is located in the waveguide at z = 0 at the upper wall. The upper 
zone (overburden) is assumed to be infinitely conducting. The depth 
dependence is such that both K and a become infinite again at a depth 
y 0 = 33km, and the conductivity of the lower zone below 33 km is also 
infinite. 

The complex relative dielectric constant may be written 

K*(z) = K 0 (l-jp)F(y 0 ,z) .(3.12) 


where F(y_,z) is the depth dependence factor 




(3.13) 


for the depth z (in km) lying between 0 and 33 km. The real relative 
dielectric constant at zero depth is K = 10 and the loss-tangent is p = 0.2, 
for the values assumed by Schwering et al (1968) at a frequency of 1 kHz. 

For these assumptions, then, the depth dependencies of the real relative 
dielectric constant and of the conductivity are 


K(z) = K 0 F(y 0 , z) (3.14) 

c(z) = 1.11 x 10" 7 F(y 0 , z) Si/m (3.15) 

respectively. If these electrical constants are independent of frequency, 





- 
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then the loss tangent at any frequency is 
p(f) = 0.2/f 

where f is the frequency in kHz. 

Schwering et al (1968) employ a z-directed vector potential to obtain 
a wave equation for this potential by using cylindrical coordinates and 
related expressions for E z , and E f . The wave equation differs from the 
usual wave equation by the addition of a term due to the ncr.-homogenecus 
depth behavior of K'(z). Going through the process of separation of variables 
for the vector potential, the authors arrive at a modal equation for each 
field component which includes "depth gain" functions. At the depth z = 0, 
the radial electric field E r = 0, and the corresponding expressions for E 2 
and H$ are simplified. 

For first three modes m = 1, 2, and 3, the attenuation coefficients at 
1 kHz are calculated to be 0.256, 1.17, and 1.99 dB/km (taken from Table II 
of Schwering et al, (1968), which lists values for the first 10 modes). 

The magnitude of the vertical electrical field strength is shown plotted 
in Figure 3.7 for 1 kHz and for the range 20 to 100 km, using the first eight 
modes. The first mode is dominant at distances in excess of 50 km; the 
effect of higher order modes is hardly evident in the range of 20 to 50 km. 

In Figure 3.8 are shown the frequency variation of the attenuation rate 
and phase velocity V i of the dominant mode (m = 1) , the latter normalized 

O 

to the velocity in free space (3 x 10 m/sec). There is resultant dispersion 
of the signal waves, but it becomes small as the frequency increases. 
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Schwering et al (1968) considered the uniform depth profile solution 

_7 

of Brown and Gangi (1963), assuming K = 10 and a = 10 Si/m for a frequency 
of 1 kHz; they obtained the attenuation for the least attenuated mode as 
0.052 dB/km. In the SPL model u.iis would correspond to the dominant (m = 1) 
mode attenuation of 0.256 dB/km. This difference emphasizes the sensitivity 
of the propagation constants to the characteristics of the assumed depth 
profile. Of importance is the effect of temperature gradient; it is zero 
for the Brown and Gangi model and about 20°C/km of depth for the SPL model. 

3.4.3 Theoretical Results of Field and Pore for More Recent Profiles 

Concerned with the possibility of lithospheric waveguide communication 
systems as a whole. Field and Dore (1973) considered the fields propagated 
via the more recent non-homogeneous profiles discussed in Chapter 2 plus 
the effects of depth attenuation of atmospheric noise and other system 
parameters such as antenna performance, receiver noise figure, and data 
rates. Frequencies considered extended from 100 Hz to 100 kHz. Using 
conceptual profiles having dry basement rock indicated more favorable com¬ 
munication transmission ranges in the 30 - 100 kHz VLF band. Presumably 
this effect was due largely to the neavier downward attenuation of atmos¬ 
pheric noise at the higher frequencies. An excellent review of their system 
analysis work was given by Bannister (1978). 

We shall discuss only their results for attenuation rates, and calcu¬ 
lations of field strength versus range and depth. 

Field and Dore (1973) used full-wave calculations for the inhomogeneous 
conductivity depth prefiles. Pectangular coordinates were employed for 
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plane TM waveguide modes propagating in the x-direction for a vertical 
electric dipole lying along the z-axis. Expressions were developed for 
the E x and Ey electric field strength components and a "voltage" normalized 
magnetic field component H’y = 120ir'Ky where Hy is the usual magnetic field 
strength. While their formulations accounted for the depth variation of 
the real relative dielectric constant, their calculations assumed a 
constant value of 10 for all depths and profiles. 

Presented in Figure 3.9 are the calculated attenuation rates for the 
modes of least attenuation as a function of frequency for the depth con¬ 
ductivity profiles of Levin (1966) shown in Figure 2.5, of Keller (1972) 
in Figure 2.9, and Housley (1973) BR profile shown in Figure 2.10. The 
least attenuated modes are TEM modes for frequencies less than about 2 kHz 
and the lowest order TM mode for higher frequencies. 

In Figure 3.9 it is seen that the attenuation depends markedly on 
both the minimum conductivity and the gradient cf the depth profile. The 
lowest attenuation occurs for the Housley BR profile, having a minimum con¬ 
ductivity of about 5 x 10' 9 Si/m. The Levin (wet) model with a minimum of 
10- 5 Si/m has the highest attenuation. While the Levin (dry) and Keller 
models have similar minimum conductivities of 10 -7 Si/m, the attenuation 
rates at 1 kHz are 0.28 dB/km for the Levin (dry) model but only 0.18 dB/km 
for the Keller profile. For the latter, the lower value is due in part to 
the sharp transition in conductivity near a 14 km depth, according to Field 
and Dore (1973). 

Note that for the uniform ideal waveguide treated by Spies and Wait 
(1971) and Wait and Spies (1971), with a conductivity of 10“ 7 Si/m for the 













waveguide, the resulting attenuations are 2 or 3 times lower than the 
values shown in Figure 3.9 for the Levin (dry) and Keller models. These 
models have more gradual "boundaries" than the sharp boundaries in the 
depth profiles of Spies and Wait. If we averaged the Keller profile over 
some range of depth, an average might be 5 or 6 x 10~ 7 Si/m; such a value 
in the Spies and Wait case would result in a higher attenuation value. 

In Figure 3.10 are presented the normalized magnitudes of the three 
field components versus depth for the Keller profile at 1 kHz. The nor¬ 
malization is such that the value of the "voltage normalized" magnetic 
field H'y is unity at the depth of minimum conductivity at 10 km. For the 
VED receiver, the magnitude of E_^ has a broad maximum near a depth of 10 km. 
Note that E x has a maximum magnitude at a 6 km depth, 14 db below the 
maximum value of E-; E x is less than E z at depths below a 4 km depth but 
E x exceeds E z at smaller depths. 

Considering the vertical field itself, the relative vertical strength 
versus frequency for several depths is shown in Figure 3.11 for the Keller 
profile. The values are normalized to that value at 10 km depth irrespec¬ 
tive of the frequency. From these curves, one notes the relative decrease 
in vertical electric field strength as the receiving antenna is moved 
toward the surface, fcr a giver, frequency. Thus, at a frequency of 1 kHz, 
the relative field strengths are 13, 32, and 44 dB below that at 10 km 
depth, for depths of 4, 2, and 1 km, respectively. At 20 kHz, the degra¬ 
dations are 16, 36, and 52 dB for those same depths, respectively. Field 
and Core have made similar calculations fcr the Housley BR model. Fcr 
depths below 10 km, the degradations in vertical electric field strength 


4 










Nounalized field strengths 


















3.26 - 


are nearly the same because the two conductivity depth profiles are the 
same. (Of course, from a system point of view, one must take into account 
the depth variation of atmospheric noise in order to obtain the resulting 
signal-to-noise ratio at the receiver. The results are not very straight¬ 
forward and one is referred to the original report (Field and Dore, 1973) 
for derivation and results. 

3.5 Other Studies - Ocean Emplaced Antennas 

3.5.1 The Work of Raytheon Co. (Grossi et al ) 

Grossi et al at Raytheon Co. (1973) had considered the fields from 

transmitting antennas emplaced in the ocean or in the lithosphere below, 

with the receiving dipole in the ocean near the ocean-lithosphere interface. 

The model employed is sketched in Figure 6.2, except that the region below 

_2 

the uniform lithosphere was assumed to have a conductivity of 10 Si/m 
instead of the infinite value assumed by Frieman and Kroll (1973). 
Frequencies of interest were in the sub-Hertz ranee. An error in the 
original expressions for the fields was corrected. There was very good 
agreement for the magnetic fields between that for a horizontal, long, 
insulated antenna on the ocean floor (which antenna was termed a "phantom 
loop" by Raytheon) and that for an insulated HED antenna when using the 
formulations of Frieman and Kroll (1973). The comparison was made at a 
frequency of 0.1 Hz and ranges of 100 and 1000 km; the lithospheric region 
had a conductivity cf 10 ^ Si/m. 

3.5.2 The Work of Frieman and Kroll 

Frieman and Kroll, 1973, determine the modes and the mode excitation 
in a flat, stratified region for both vertical and horizontal current 





distributions. The region of interest lies between perfectly conducting 


plane boundaries at z = -d^ and z = h where d^ is the depth of the litho¬ 
sphere and h is the height of the ionosphere, both measured from the ocean 

floor.* The depth of the ocean is d w . The parameters used are j w = 4 Si/m, 

-6 

o ^ = 10 Si/m. With frequencies in the range between f = 0.1 and 10 Hz, 
the propagation constants are, for the lithosphere Y ? = (jwuc^) 15 , for the 

L 2 1* 

sea water, Y w = (jwua w ) and for the air, Y a = (-<d ue ) . By confining 

their attention to modes that propagate to distances large compared to both 
dj and h, Frieman and Kroll require only two modes. These are the litho¬ 
spheric mode which is large in the lithosphere and the air mode which is 
large in the air. Solutions for the magnetic field of these modes are ob¬ 
tained with both vertical and horizontal electric dipoles located on the 
floor of the ocean. Their results indicate "that at large distances the 

HED is more effective than the VED bv a factor (c /a.)* 5 = 2x10^", They go 

w ** 

on to show that the contribution from the terms that involve the depth d^, 
of the ocean are essentially the same as for an ocean of infinite depth. 
Finally, in the range where the thickness d^ of the lithosphere is effec¬ 
tively infinite, the magnetic field reduces to that obtained by Banos, 1966, 
or Kraichman, 1970, for the problem of two half-spaces. This is given by 
the near-field formula 


He. 


_ cos - 

2*Y.. „ 3 


In a specific example, the magnetic field is 
insulated antenna 10 km long and 4 cm in diameter 


calculated for a terminated 
lying on the ocean floor. 


See Figure 6.2 for geometry employed 





4. The Two-Half-Space Model and Lateral-Wave Transmission 
4.1 Introduction 

In their modal analysis of electromagnetic wave propagation in a litho¬ 
sphere composed of granitic rock (j= 10 ° Si/m) and bounded above by sea 
water (o ■ 4 Si/m) and below by the mantle (o = « ) Frieman and Kroll [1973] 
concluded that with horizontal-electric-dipole excitation and reception at 
points on the ocean floor, the received field was very little affected by 
the thickness of the lithosphere (d„ "v 10 km) or by the depth of the water 

(d ^ 5 km) in the frequency range from 0.1 to 10 Hz. This means that the 

w 

electric field at the point of reception is governed by a single dominant 
mode which corresponds to the solution for transmission along the boundary 
between two half-spaces, the ocean and the lithosphere. This conclusion 
is confirmed by the careful interpretation by Bubenik and Fraser-Smith [1978] 
of their analysis of the field of a vertical magnetic dipole located in 
sea water midway between the air above and the sea floor (lithosphere) below. 

Bubenik and Fraser-Smith compared the contributions to the total field by 
direct transmission through the sea as if this were infinite in all directions 
with those by air-surface and sea-floor lateral-wave modes. They concluded 
that for a dipole near the air surface, the air-surface mode dominates; for 
a dipole near the sea floor, the sea-floor mode dominates. Since a vertical 
magnetic dipole is equivalent to four horizontal electric dipoles arranged 
in a square, its field is a superposition of the fields of the four dipoles; , 
each like the field of a single horizontal electric dipole. The superimposed fields 
of the very closely spaced dipoles lead to a large reduction in the resultant field 
due to cancellation but this in no way alters the nature of the field with 
respect to its modal structure. This means that if the sea-floor lateral-wave 
mode dominates for a vertical magnetic dipole in sea water near the sea floor, 
the same is true for a horizontal electric dipole at the same location. 
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Extensive numerical calculations by King and Sandler [1977] based 

on the exact general theory of Banos (1966] have shown that transmission 

between horizontal electric dipoles in a dissipative half-space near its 

boundary with air is predominantly by the so-called lateral waves. This is 

true for a dissipative half-space that has the properties of sea water, 

9 

lake water, or dry earth over a frequency range from 1 to 10 Hz. For 
radial distances from o = 1 m to 100 km and more, the dominant component 
of the electric field at a depth z in the medium due to a transmitting dipole 


at the depth d is the radial one, i.e., E . While there are short ranges 

p 

where the component E, is somewhat greater than E , this difference is 

9 P 

never large. On the other hand, there are large ranges including especially 


the greater distances where E greatly exceeds E . It follows that for 

P $ 

practical transmission over large distances E is the only satisfactory 

P 

component. It has also been shown by direct comparison of the numerical 


calculations from the exact integrals with a set of approximate formulas 
due to Banos [1966] that these latter are roughly correct at least in order 


of magnitude far beyond the ranges of validity specified by Banos. Specifically, 
the graphs of the electric field as a function of distance at a given frequency 
as calculated from 3anos' near-, intermediate-, and far-field formulas can be 


extended to points of intersection. The sequential combination of the three 

graphs provides a continuous curve of IE 1 against distance from the source o 

‘ P 

(but with discontinuities in slope at the intersections of the near- and inter- 
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mediate-field graphs and of the intermediate- and far-field graphs) that 
follows fairly closely the smooth graph computed from the exact integrals. 
This general agreement is very significant because it justifies the use of 
Banos' simple formulas to gain an understanding of the physical nature of 
the waves propagating from the source to the receiver and of the significance 
of the numerous parameters and variables that combine to determine the field 
at the point of observation. 


4.2 Banos' Approximate Formulas and Their Interpretation 

The coordinates and parameters used in the representation of the field 
in a region consisting of a half-space of sea water (Region 1, z > 0) and a 
half-space of the lithosphere (Region 2, z < 0) are shown in Fig. 4.1 . The 
wave number in sea water with the time dependence e lwt is* 

2 1/2 1/2 
k l “ S 1 + ict l = (w u 0 £ l + iu Vl } 4 (i “ li 0 a l ) (4.1) 


The wave number for the lithosphere is 


2 1/2 

k., = 3^ + la ^ = U> y o £ 2 + 


(4.2) 


As indicated, the sea water behaves like a good conductor over all useful 

frequencies. The same is true for the lithosphere at low frequencies. At 

high frequencies the lithosphere may behave like a dielectric so that the 

more general formula for the wave number is retained. 

Bafios’ approximate formulas for the electric field of a horizontal electric 

dipole at the location shown in Fig. 4.1a are given by (7.10), (7.32) and (7.45) 

in his book [Banos, 1966]. With small changes in the notation and the substi- 
2 

tution of k“/iuiUQ for o^, the formulas for the radial component of the electric 
field are: 


Wear Field: 


'lo 



cos ^ 
3' 

0 


i[k 2 P + k^(d+z)] 
e 


(4.3) 


Intermediate Field: 


*The form on the right in (4.1) assumes a./uie^ >> 1 for sea water. 
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4.1a Two-half-space Model for Lateral-Wave Transmission 
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lb Schematic Diagram of a Ray in Lateral-Wave Transmission and 
its relation to the Mathematical Representation: 
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io)U r .k^ cos 5 k. . i [k ? p + k.(d+z)] 

E. A-7-=-- [1 + t — (iik p/2) 1/Z ] e " 1 (4.4) 

?n\c~ 0 k. 


Far Field: 


- • MU 0 cos > Hk 2 p + k x (d+z)] 

L lp ’ 2wk. 2 e 

2 0 
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Tae formulas for the other components are similar, but since is the only 
generally useful one, it is sufficient to consider it. The ranges in which 
each of the three formulas is a good approximation have been specified by 
Banos. For sea water they are extensive, but even for this they do not actually 
join. Numerical calculations from the exact integrals have shown the three 
formulas may be used successively as the radial distance is increased with 
the transition from the one to the other occurring where the magnitudes of the 
field are equal. At and near these transitions from the one formula to the 
other the approximations are quantitatively inadequate but qualitatively 
correct. 

Tlie lateral-wave nature of the field is evident in all three formulas. 

Note that they all have the same exponential term which is readily interpretted. 
It indicates a traveling wave originating at the source dipole, (p=0, p=0,z=d), 
traveling downward in the sea water to the boundary surface , then traveling 
radially outward in the lithosphere a distance p , and finally traveling upward 
from the surface to the receiver at the point of observation at P(p,5,z). This 
involves travel in the sea water (Region 1) over a total distance z + d and 
travel in the lithosphere (Region 2) over the distance o rFig. 4,1b). In the sea 
the attenuation and phase constants are, respectively, a 1 and 3^; in the litho¬ 
sphere they are and If the real and imaginary parts of the exponents 

are separated, the common exponential term has the expanded form 


-[c^P + a^(d+z)] i[3 0 p + 3^(d+z) 

e e 


(4.6) 




Clearly there is attenuation in both regions but since the conductivity of the 
lithosphere (Region 2) is much smaller than that of the salt water (Region 1) 
it follows that << a^. For a given range p a n important quantity in deter¬ 
mining the amplitude is the total vertical distance z + d. The closer the trans¬ 
mitting and receiving antennas are to the boundary, the larger will be the ampli¬ 
tude of the field. However, if a^(d+z) < 1, g^^W+z) ^ ^ even for large 

values of d + z. This can be accomplished by choosing a low enough frequency so 
that is small even though is as large as 4 S/m. Insofar as the exponential 
term is concerned, small distances z and d are desirable as well as a frequency 
that is not too high. For a specified range the contribution to the attenuation 

“Ct 0 

by the factor e 2 can be reduced only by making ct 0 smaller. With a given conduc- 

tivitv of the lithosphere a 0 can be decreased only by lowering the frequency. 

If the exponential term does not dominate,a significant field can be 

maintained at great distances. In order to understand its dependence on the 

distance and frequency it is necessary to examine the other factors that 

determine the amplitude. The radial distance appears in all three formulas 

(4.3) - (4.5) but in a different manner in each. In the near-field formula 

3 

the amplitude decreases as l/o , a prohibitively rapid rate of decrease. 

In the far-fiela formula the rate of decrease is as l/o~ , a great improvement 

over the near field. Finally, in the intermediate range the rate of decrease 

with radial distance varies between l/o and l/»p - which is by far the most 

favorable. The comparable dimensionless quantities for the three ranges are: 

-3 2 3 2 

(k^o) , k“/k^o , and l/k^k ? p in the order near field, intermediate field, and 

far field. With a proper choice of frequency for the parameters of the media 
involved, the largest magnitude of the radial electric field should be obtained 
when o is in the intermediate zone near its boundary with the far zone. This 
clearly excludes very low frequencies for which the near-field formula is appli¬ 
cable for the relevant radial range. 

A better understanding of the complete picture including the exponential 
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and inverse distance dependences can be obtained from Fig. 4.2 which displays 

contours of constant IE (in d3 referred to 1 volt per meter) as a function of 

P 

the radial distance and the frequency for a horizontal electric dipole at a 
depth a = 15 cm in water. The point of observation is at the depth z = 30 cm in 
the water. The half-space bounding the water is air filled. In the previous 
notation, the water is Region 1, the air is Region 2. In this case 012 = 0 and 
$2 = 2 tt /\2 where X 2 is the wavelength in air. The real effective permittivity of 
the water is e , = 30; the real effective conductivity- in Siemens per meter is 
a ^ = 0.004, 0.04, 0.4 and 4 in the diagrams from left to right. The diagram 
on the left applies when the antennas are in lake water of high purity, the dia¬ 
gram second from the left is for polluted lake water, the diagram on the extreme 
right for sea water, The contours were obtained from numerical computations 
of the electric field using the exact integrals. 

In the diagram on the left the conductivity is sufficiently low so that 

the exponential attenuation, e is ll0C dominant at any frequency with 

z + d = 45 cm. In terms of the approximate formulas of Banos, the right-hand 

third of the diagram is for low frequencies and the near-field formula with 

,3 

its l/o amplitude factor, the left-hand third is for high frequencies and the 
far-field formula with its 1/p^ amplitude factor. The transition between them 
at intermediate frequencies is for the intermediate-field formula. Note that at 
any fixed distance, for example c = 8 km, is very much greater at the 

high frequencies than at the low ones. Specifically, at the high frequencies 
E, = -120 dB, in the low-frequency range E. = -200 dB. In the second 

-*-C lo 

diagram from the left the conductivity of the water is ten times greater but 
still small enough to make the exponential attenuation unimportant. The contours 
of constant electric field are quite similar to those on the left with only a 
small decrease in amplitude in the high-frequency range, a somewhat larger one 
in the low-frequency range. In the third diagram from the left the conductivity 
is multiplied by another factor ten and this makes it large enough so that the 
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exponential attenuation begins to play a significant role, but only at the 
nigner frequencies. The rapid decrease in amplitude at frequencies above 10^ Hz 
is due to the exponential factor e 3 -^ 2+c ^^ with z+d = 45 cm. Note, however, 
that even with the exponential attenuation significant, the field in the high 
frequency range is still much greater at any fixed distance than in the low 
frequency range. The optimum frequency is clearly in the vicinity of 10 / Hz 
where the contours show a maximum distance for any specified value of E j. 

In Che graph on the right for sea water the exponential attenuation has 
become the major factor at the higher frequencies. The radial electric 

field decreases to extremely small values at all distances when the frequency 
g 

is increased to 10 Hz or higher. However, at any fixed distance the field 
has a maximum near f = 10 / Hz and this maximum is significantly greater than 
the field at the same distance at any of the lower frequencies. Thus, when 
z + d = 45 cm the field at any radial distance p is greater at a frequency 
between 5 x 10^ and 5 x 10 7 than at any frequency below f = 10^ Hz. 

If the total distance of travel in the sea water is increased from z+d 
= 45cm the maximum field at any fixed radial distance moves toward lower 
frequencies. When this reaches about 10 3 Hz the exponential attenuation will 
effectively eliminate the maxima in the high frequency range and the field 
in the low-frequency range will be larger than that anywhere in the high-fre¬ 
quency range. The choice of frequency for lateral wave transmission is seen 

mmr X (z+d ) 

to depend critically on the quantity i^(z + d) and the exponential e 1 
When the region of low conductivity is not air with a, = 0 but the lithosphere 

-8 -4 

with o 0 = 10 to 10 .the additional exponential attenuation by the factor 
e must also be considered. This requires a complete recalculation. 

However, the basic properties of wave propagation with horizontal-dipole 
excitation near an interface are best understood from the simpler case with 
a perfect dielectric like air as the half-space with small conductivity. 

3efore proceeding to a study of lateral ’.rave propagation along the 
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surface of two half-spaces with the properties of sea water and the lithosphere, 
it is well to point out that the simple picture illustrated in Fig. 4.1b of 
lateral waves traveling vertically downward from the source dipole a distance 
d to the boundary surface, along the surface a radial distance p to the point 
below the receiver, and then vertically upward a distance z to the receiving 
dipole, is an idealization that is an excellent approximation with a highly 
conducting medium like sea water. With poorer conductors the optimum angle 
downward from the source to the surface is the critical angle at which the 
incident wave is neither transmitted nor totally reflected but travels along 
the surface. With sea water this angle is very close to 90°; for media with 
lower conductivities and permittivities the angle can be considerably smaller 
than 90°. A discussion and determination of this angle is given by Staiman 
and Tamir [1966]. For a maximum transmitted signal the dipole should not 
be parallel to the interface but perpendicular to the direction of the critical 
angle. With sea water the difference between the two orientations is negligible. 
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4.3 Lateral-Wave Transmission Along the Lithosphere-Ocean Boundary 

Lateral-wave transmission along the boundary between air and water involves 
exponential attenuation only along the two vertical parts of the path, both in 
the salt water. This has the form e -0 * 1 ( z+ ^) w h ere j the depth of the horizon¬ 
tal dipole in the sea water and z is the depth of the point of observation. When 
the transmission is along the boundary in air it is loss free since e^ 20 = e*""-" 
with k,, = + ia 0 and = 0. When Che half-space of air is replaced by a naif- 

space of rock with the real effective permittivity e g? and the real effective 

conductivitv a the wave number k„ = 8„ + ict„ is complex so that e^ 20 = 
e2 2 2 2 

e i32P • Thus there is exponential attenuation along the main path of 
propagation in the lithosphere in addition to the decrease in amplitude governed 
by the complicated functional dependence on the inverse powers of the radial 
distance p. 

Computations have been made from the exact integrals for lateral-wave 
transmission from a horizontal electric dipole with unit moment, Id£ = 1. 

The dipole is at the depth d in the sea water and the field is calculated at 
the radial distance p at the depth z = d. (Note that in the computations for 
air and sea water, z = 2d = 0.3m) . The material parameters used in the cal¬ 
culations of the electric field are e , = 80e„ and o , = 4 Si/m for Region 1, 

el 0 el 

—8 

the sea. For the lithosphere, Region 2, t 7 = 16 eq and ranges from 4 x 10 

Si/ra to 4 x 10 Si/m. Initially z = d = 0.15m . Computations of !E ' (which 

1 0 

is the dominant component in lateral-wave transmission) have been made over 

9 

radial distances from p= 0.1 to p = 100km at frequencies from 10 to 10 Hz. 

Graphs showing 1 E ! as a function of the radial distance with the frequency as 

P 

—8 —7 

the parameter are shown in Fig. 4.3a for o = 4 x 10 and 4 x 10 Si/m and 
in Fig. 4.3b for a , * 4 x 10 ^ 4 x 10 \ 4 x 10 \ and 4 x 10 ^ Si/m. From 

these graphs the contours of constant |E ! shown in Figs. 4.4a - 4.4c have been 

P 

constructed. These show E ! as a function of both the radial distance and 

P 

the frequency in the manner carried out in the comparable Fig. 4.2 for air and 


water. 
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The set of contours on the left in Fig. 4.4a is for sea water and a litho- 

—8 

sphere with the low conductivity a= 4 x 10 Si/m. It resembles the correspon¬ 
ding set of contours on the extreme right in Fig. 4.2 for sea water and air with 
the conductivity c? e? = 0. In both sets of graphs the exponential attenuation 
e -M dominates the high-frequency quarter on the left. Here the graphs show 

a rapidly decreasing radial range for any given value of j E | as the frequency 

P 

is increased. The rate of decrease for sea water and air is even more rapid than 

that for sea water and the lithosphere because (z+d) = 0.45m with the former, 0.3m 

with the latter. On the other hand, attenuation by e 1 is insignificant in 

the low-frequency third of the contours on the right in each diagram. This is 

because at low frequencies <*2 is very small. For air, a 2 = 0, e -020 = 1, and 

the contours for fixed values of are virtually horizontal in the low-frequency 

third of the right-hand set of contours in Fig. 4.2. This indicates that the 
radial range for a specified field is almost independent of the frequency in the 
low-frequency range. For the lithosphere, U 2 ^ 0 and e a “° varies significantly 

with both the frequency and the radial distance. As a consequence, the contours 
for given values of E at low frequencies range from horizontal lines when c^ 
is small to steeply increasing graphs when a ~ is large. The increase in radial 
range as the frequency is reduced and the conductivity increased is clearly shown 
in the sequence of graphs in Figs. 4.4b and 4.4c. 

Perhaps the most significant characteristic of the contours of constant 

, E ; for sea water and air (on the extreme right in Fig. 4.2) and for sea water 

P 

and the lithosphere (Figs. 4.4a - 4.4c) is the relative maximum in the radial 

range which occurs in all diagrams at frequencies between f = 1 and 10 MHz so 

_3 

long as a ,,< 4 x 10 Si/m. Particularly notworthy is the fact that the rela¬ 
tive maximum value of p for a given E_ is much greater than the corresponding 

distance for the same value of E : at very low frequencies when a^ * 4 x 10 

“8 

Si/m. This can be seen from Figs. 4.4a and 4.^b for both a^ = 4 x 10 

t , « 4 x 10 Si/m. For example, when p , 5 4 x 10 Si/m, the value 


and 






E ; * -175 dB occurs at c 3 40 km when the frequency is between 10^ and 10^ Hz 

but o < 0.3 km for all frequencies from 10 to 10"* Hz. Similarly [E [ = -225 dB 

occurs at o near 500 km when the frequency is between 10^ and 10^ Hz but at 

4 -7 

2 km for all frequencies between 10 and 10 Hz. When a „ = 4 x 10 Si/m the 

e2 
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graphs are very similar to those with a 0 = 4 x 10 Si/m with a reduction to 
p ■ 100 km in the radial range at the maximum between f = lO* 3 and 10 ' Hz and 
virtually no change in the short range at low frequencies. However, when o g0 
is as large as 4 x 10 b Si/m the maxima in the radial range near f = 10^ Hz are 

substantially depressed but the radial distance for the same values of Je ' at 

P 

very low frequencies remains smaller so long as p < 75 km. Thus, |E I = -225 

o 

dB occurs at p * 33 km at the maximum with the frequency between 10° and 10 ‘ Hz 

3 

and drops to p = 2 km at f = 10 to 10 Hz. For radial distances greater than 

about 100 km, a given value of |E | is found at substantially greater distances 

P 

at the very low frequencies than at the relative maximum between f = 10° and lo" 

Hz. This is a consequence of the radial exponential attenuation factor e 

P 7 

which is small enough at frequencies as high as 10 and 10 Hz to depress the 
relative maximum more and more as p is increased. This effect is dramatically 

shown in the four sets of graphs in Figs. 4.4b and 4.4c in which the conductivity 

”6 * 3 

of the lithosphere increases from 4 x 10 on the left in Fig. 4.4b to 4 x 10 

n 

Si/m on the right in Fig. 4.4c. It is seen that E \ at the maximum near f = 10 

p 

Hz is reduced more and more while |E_J at very low frequencies and sufficiently 
great radial distances is increased. 

The contour diagrams in Figs. 4.4a,b lead to the following general conclu¬ 
sions regarding lateral-wave transmission along a plane boundary between sea water 
and the lithosphere when the source is a horizontal electric dipole in the sea 
at the small distance d = 15 cm from the bounding surface and the point of obser¬ 
vation is at the same distance z = d from the surface at a radial distance o : 
a) If the conductivity of the lithosphere is greater than * 4 x 10 ° Si/m, 
transmission over distances greater than 50 km experiences the lowest possible 
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attenuation at frequencies below 1 kHz. This is true for any values of z = d 
>15 cm. 

b) If the conductivity of the lithosphere is significantly smaller than 4 x 10 ^ 
Si/m, there is an optimum frequency between 1 and 10 MHz at which the radial 
range has a relative maximum that may be substantially greater than the radial 
distance for any value of |E | at other frequencies including the very low. 
Furthermore, an effective directional array with significant gain can be con¬ 
structed for both the transmitter and the receiver for use at frequencies in 
the 1 to 10 MHz range. 

c) The conclusion under b) assumes the transmitting and receiving antennas to 

be in the sea water but very near the boundary surface with the lithosphere 

(d = z = 0.15m). It is to be expected that as the distances d and z from the 

antennas to the surface is increased, the exponential attenuation e - 'M( z+t ^ 

will become dominant at lower and lower frequencies. This must simultaneusly 

depress the relative maximum which occurs between f = 1 and 10 MHz when (z+d) = 

30 cm and move it to lower frequencies. In order to investigate this effect 

quantitatively.computations were made of 'with z = d = 1.5 m over the same 

range of conductivities for the lithosphere as used in Figs. 4.3a and 4.3b. 

These results are shown in Figs. 4.5a and 4.5b. From them the contour diagrams 

in Figs. 4.t>a - 4.6c were constructed. They correspond to the diagrams in 

Figs. 4.4a - 4.4c. It is seen that when a , < 4 x 10 3 Si/m the relative 

maxima in the radial range have moved from near f = 10^ Hz wit 1 ' z * d = 0.15m 

to near f = 10 3 Hz with z = d = 1.5m and that they are substantially smaller. 

—8 

For example, with a = 4 x 10 Si/m the maximum range for |E » -225 dB is 
p = 10 km with z = d = 1.5m instead of p v 500 km with z = d = 0.15m. There is 
no corresponding decrease in the range at low frequencies since there the expo¬ 
nential attenuation is small. It follows that unless both the transmitting and 
receiving antennas are very close to the ocean-lithosphere boundary, the greatest 
radial range is achieved at low frequencies. 
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d) There is a potentially useful low-attenuation window in the frequency range 
between 0.1 and 10 MHz. It exists only when the conductivity of the lithosphere 
is sufficiently low, \ x 10 ° Si/m. Furthermore, it is available only 

when the transmitting and receiving antennas are very close to the interface 
between the ocean and the lithosphere. When (z-M) = 0.3m the maximum radial 

range occurs when the frequency is near 10^ Hz. As the distance from the 
boundary is increased the frequency for the maximum range decreases. It is near 
f = 10^ Hz when (z+d) = 3 m. The magnitude of the maximum range also decreases 
as (z+d) is increased. 
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5. LABORATORY MODELLING AND ELECTRODYNAMIC SIMILITUDE 

5.1 Introduction 

Laboratory modelling has been used to simulate or demonstrate the 
behavior of large scale electromagnetic systems by scaling dimensions and 
frequencies while keeping the same shape of the model and the system it is 
to simulate. For example, a full scale linear antenna of length L g operating 
in air at a frequency f can be simulated by measurements on a similar shaped 
antenna of length = L e /n provided the model is operated at a frequency 
f m = n f e in air. For more general lossy media, a knowledge of the prin¬ 
ciples of electrodynamic similitude and limitations is needed for the use 
of models and interpretation of data. We examine relations for application 
to some propagation problems. 

Millington (1948, 1949a,b) discovered the "recovery effect" on the 
amplitude of LF and VLF ground waves as they traversed a land-sea boundary. 

A simple model was used to demonstrate the effect, wood to represent the 
land and copper sheeting to represent the sea. A small microwave source was 
moved on a small cart away from the source on wood onto the metal whence the 
signal amplitude increased markedly at the "coast-line," (recovery) and then 
attenuated slowly along the copper "sea". Quasi-empirical relations were 
developed in his publications.* Pressey (1956) developed quasi-empirical 
relations for the phase shift. The matter was studied theoretically and 
experimentally by King and Maley (1966) and King, Maley, and Wait (1967). 

*A group of us in February, 1950, witnessed Millington's demonstration 


at Marconi's arranged through the courtesy of the late Sir Robert Watson Watt. 
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The model was a tank filled with water to represent land and aluminum sheets 
to represent the ocean; data were obtained at 4 GHz. While scaling was not 
exact, the measured electric field strength was reasonably well in accord 
with prediction. 

Model studies of propagation in an idealized crustal waveguide were 
made experimentally and theoretically by Brown and Gangi (1962, 1963) and 
Gangi (1966). We shall discuss their results, arranged to demonstrate 
similitude, a bit later, for a uniform slightly lossy guide represented by 
saline water. 

Iizuka (1967, 1968, 1969) studied the impedance of an antenna and the 
field variations of waves propagated in a lithospheric waveguide for a guide 
having variable conductivity with depth. A gel of agar-agar with diffusion 
of sodium chloride was arranged to provide the desired profile. Similitude 
in size, frequency, etc., was accomplished and results will be discussed a 
bit later, also. Others have developed the similitude relations. 

5.2 Principles of electrodyr.amic similitude 

The principles of electrodynamic similitude were developed by Stratton 
(1941) and more completely delineated by Sinclair (1948). Brown and Gangi 
(1962, 1963) and Gangi (1966) gave the resulting equations, based on 
Stratton (1941) relating the scaling of frequency, sizes, and electrical 
parameters of the model system in terms of those quantities of the original 
system. Keller and Fritschkneckt (1966) gave a brief discussion for 
scaling media of large loss tangent. Iizuka (1968) developed the relations 
while employing generalized orthogonal coordinates for the identities of 
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Maxwell's equations in the original and model systems. Fritschkneckt (1971), 
basing his work in part on that of Sinclair (1948), used rectangular coordin¬ 
ates as did Sinclair (1948). 

Since we are dealing with systems immersed in the earth , such as the 
lithosphere, we use the subscript "e" on related quantities; the subscript 
"m" will denote those quantities for the model . 

The original and model media are characterized as isotropic, homogeneous, 

or inhomogeneous, with electrical constants of permittivity c = e Q e r (where 

—9 

e Q is approximately 10 /36n F/m and e r is the relative dielectric constant, 

a permability y (which we will ultimately assume to be that of free space 
u D = 4 it x 10 -7 H/m) and a conductivity a in Si/m. The unit of length is r. 

In order that Maxwell's equations remain invariant from the original to model 
system, we shall scale the electric field strength E, the magnetic field 
strength H, and time. Thus, for the fields: 
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wave frequency and w the angular frequency. The equations (5.1 through 5.3) 
are cast in the form of Iizuka (1968)*; as mentioned by Fritschkneckt (1971), 
they are termed by Sinclair (1948) as "absolute" models, for which q E and q H 
are known separately and not merely their ratio. If a specific value is not 
assigned to q E or q^, the model is termed "geometric" by Sinclair (1948). 

Multiplying equation (5.3) by (5.4), one obtains: 


r_ u)_ vuZcZ 
m in m m 


w e ^"e £ e 


(5.6.a) 


k_'r = k 'r 
mm e e 


r /X ' = r /X ' 
nr m e e 


(5.6.b) 

(5.6.c) 


where k ' = wj/ne = 2ir/X ' and k ' * w v/yTeT = 2tt/X„ ' ; i.e . , the k ' 
m irr m m m e eee e — — m 

and k ' are the wave numbers in loss-less dielectrics.** 
e 


From (5.4) and (5.5): 


(5.7.a) 


or 


U) € (jl) € 

mm e e 


60o X 


2 = p = 
r e 


60*7 X 

e o 


(5.7.b) 


rm re 

where d is also known as the loss-tangent tan 5 and X„ is the free 
m,e m,e° 

space wavelength (meters). 


* i.e. , equations (IV-13) through (TV-15) of Iizuka (1968). 

♦♦Equations (5.6) through (5.9) correspond to Iizuka's (1968) equations 


(TV-16) through (IV-19) 
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4 

A* 


Again, from equation (5.3) 


— r a u = r w u 
in m m e e e 


(5.8) 


^E 


which can be changed,* by using (5.6), to show: 


? E fee l^m 


q„ G 1 U 1 £ 
h H m 


(5.9) 


rm e 


This relation is used when showing the scaled impedance relation of original 
and modelled antennas. 

lizuka (1968) preferred as convenient to express the similitude condi¬ 
tions as follows: equation (5.6) states the identity of electrical sizes in 
terms of the loss-less wavelength between original ("e") and model ("m") 
systems, and equation (5.7) states the identity of loss-tangents cf the two 
systems. 

Brown and Gangi (1962, 1963) and Gangi (1966) gave relations similar to 
those above, viz.: 


f u e r -f u e r 
m m m m e e e e 


f m u m °m V = f e u e e e V 


(5.10) 


(5.11) 


If one takes the ratio of (5.11) to (5.10), the identity of loss-tangents 
in equations (5.7) result. 


•Equations (5.6) through (5.9) correspond to lizuka's (1968) equations 
(IV-16) through (IV-19). 
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With the assumption that U m = U g = U Q , they go to show: 



(5.12) 


£m 

f d 



(5.13) 


where D e is the depth of the original lithospheric waveguide, and D m that 
for the model. While they do not give expressions similar (5.9), Brown and 
Gangi in their published article (1963) do note that scaling keeps the 
number of wavelengths between scaled distances constant, which is equivalent 
to equation (5.6.c), i.e ., loss-less wavelengths. 


In practice, there are severe restrictions on using these relations in 
a simulator model. Paraphrasing Iizuka (1968): 

(a) the values of e r and loss tangent for the model are limited 

(b) the dimension scale factor has to be such that the model is not 
too small or too bulky 

(c) ranges of e r and loss-tangent by mixing materials is limited 

(d) changing loss-tangent or e r is complex, it being possible to 
minimize one while changing the other 

(e) materials must not be noxious, explosive or too costly. 


3 laboratory lithospheric models and similitude 

3.1 Idealized waveguide with electrical constants uniform, with depth 

Brown and Gangi (1962, 1963) and Gangi (1966) made laboratory measure¬ 
ments with a simulated waveguide, the simulation representing an "average" 







set of conductivity values deduced from their proposed profiles (See 
Chapter 2). 

The model of the lithosphere to be simulated is shown in the Figure 5.1. 
In their measurements, the earth model was assumed to have a depth D e = 35 km, 
with e r 2 = e re = 4. Since their conductivity depth model (Chapter 2) showed 
a conductivity $2 ~ might vary with frequency, appropriate scaling 

ratios were employed. 


An aluminum (or carbon) plate represented the overburden or upper crust 

and an aluminum plate was used for the lower crustal boundary (Moho). The 

crustal waveguide was a salt water solution whose conductivity could be 

5 

varied. A frequency scale factor of f m /f e = 10 was chosen for simplicity 

in the models. The value of c m for water was 80. Since waveguide distances 

of 1000 km were desired to be simulated, and a conveniently-sized tank was 

desired for the laboratory, a distance (or depth D m /D e ) scale factor works 
5 

out to be 4.47 x 10 ; thus if D e were 35 km, a tank depth of 7.82 cm or 
3.08 inches results. Distances were plotted in units of H g or 

< H e " D e' % = D m>• 



In Figure 5.2 is shown measurements on the model at 300 MHz, in the 

tank 3.08 inches in depth, of the vertical (relative) electric field strength 

over a range of 2H to 17H (when H is H m or H g , the conductivity a m being 

0.07 Si/m. The loss tangent p m = .0525 which is small; the corresponding 

/, -6 

plane wave attenuation is then = 60iro m A /e^ =* 3.30 x 10 n/m = 0.28 
dB/km, These values of loss-tangent for the earth p e and a g work out to be 
the same, since they are scaled. 
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In comparison with theory, Figure 5.2 shows bracketing calculations 

as dots or crosses; the dots correspond to 0 m = .05 Si/m for which the 

loss tangent p_ is 0.0375 for f = 300 MHz just as it is for the scaled 
in m 

-8 

value of a g = 2.5 x 10 Si/m at a frequency of 3KHz. The crosses are 

for calculations for a = 0.10 Si/m for which the loss-tangent p = 0.075 

m ^ ^m 

—8 

at 300 MHz, corresponding to the earth values = 5 x 10 Si/m at 3 KHz. 

The calculations bracket well the measured curve. 

Brown and Gangi (1962, 1963) and Gangi (1966) show the effects on 
field distributions when simple inhomogeneities were introduced. These 
included (a) a step barrier in the lower boundary to correspond to the 
shallower depth of the Moho under the sea at the continent-sea boundary; 

(b) a bar in the upper boundary to represent a geologic fault zone, and 

(c) a cylinder to correspond to a lens or intrusion of mountains into the 
granite. While idealized, the results showed potential interference pat¬ 
terns compared with those when no discontinuity was present. (See the 
original papers for detailed results). It should be noted that the ocean 
crust may differ markedly in electrical properties from that of the continen¬ 
tal crust (Cox, 1971, 1978). 

5.3.2 Lithospheric waveguide with constants varying with depth (inhomogeneous ) 

An inhomogeneous lithospheric waveguide, i.e ., a waveguide whose elec¬ 
trical parameters vary with depth has been investigated theoretically in 
the past seven or eight years but relatively little has been reported on 
simulation experiments. While the depth variations may be represented by 
thin horizontal slabs, it is a bit more satisfying to simulate the smooth 
(sic) depth profiles by suitable materials whose properties may be smoothly 
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changed. About the only extended model investigations reported are those 
of Iizuka (1967, 1968, 1969, 1971, 1972). I have found no theoretical 
comparisons with his measurements. His results extend earlier ones on 
antenna impedance with stratified films of Iizuka and King (1962 a,b). 

The simulator chamber used an agar-agar solution (gel) into which a 
solution of sodium chloride (NaCl) could diffuse to change the concentration. 
Two chambers were built, one for antenna impedance measurements and a second 
larger one for field strength measurements. The technique of fabricating 
inhomogeneous media was described in Iizuka*s 1967 paper; a summary of the 
work on impedance and field strength appears in his 1972 paper. 

The variation of loss tangent and relative dielectric constant as a 
function of frequency are shown in Figure 5.3 for various media such as 
agar-agar, agar-agar with ethyl alcohol, tap water, 0.1 mol NaCl solution 
and ethyl alcohol. For agar-agar NaCl affects loss tangent; ethyl alcohol 
affects e r . 

Agar-agar has a high dielectric constant and for simulation use it 
can reduce the dimensions of the simulation chamber. For various reasons 
(See Iizuka 1968, Chapter III), a modelling frequency of 114 MHz was chosen. 
The properties of agar-agar with concentration of NaCl at 114 MHz are shown 
in Figure 5.4. On the Figure equations for least squares fit are given for 
e r and tan 6, as a function of the concentration C of the NaCl. The agar- 
agar solution was 1.3% agar-agar and 98.7 tap water. The concentration of 
NaCl was the weight in grams in 1000 cu cm of agar-agar. If this concen¬ 
tration C is less than 3, then the approximations simplify to: 




Oielect.ric constant and Loss-tan er.t vs. frenuency for various neLin,especially a ,T ar-a rr ar 

(fro-. Iiy.uka, I 06 rt ) 
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electric constant and loss •tangent of a'ar-a-’ar as a function of concentration 
of sodiuo chloride (Nall) at 114 MHz. (Fro - Iizuk?., 1468 ) 
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e r - 78.30 + 0.073 C 


p = tan 6 =* 0.1569 + 0.4633 C 


(5.14) 


(5.15) 


for 114 MHz. 


A choice of f m can be made by consideration of Figure 5.3, and equation 
(5.6) or (5.13). From Figure 5.3. e r of agar-agar is near 78 for frequencies 
less than 4000 MHz, then decreases at higher frequencies. The loss-tangent 
starts to increase near 4000 MHz. This limits the agar-agar choice of n f m , 
n being the dimension scale factor r e /r m , and the original full scale fre¬ 
quency f . Thus from equation (5.13): 


f n/ f e " n vz re /z r 


(5.16) 


The relative dielectric constants of the agar-agar e m vs f m are known from 


Fiqure 5.3. In Figure 5.5, the model simulation frequency f m is plotted 
vertically against the product of the full scale frequency f e and the dimen¬ 
sional scale factor n for various values of relative dielectric constant E re 
of the full scale system. The lowest loss-tangent of agar-agar is that get 
with no NaCl. It is plotted on the left of Figure 5.5 vs the model frequency 
^m* These curves can determine model parameters. If the lithosphere wave¬ 
guide is to be operated at f e =5KHz, with e re * 4, and if the distance factor 
5 

is n * 10 , then n f e = 500 (if f e is in Mhz). The model frequency then is 
f m = 114 MHz at which the minimum loss-tangent is 0.13. 


Various profiles of the lithosphere itself were deduced from values pub¬ 
lished up to about 1965 and combined into one plot from Iizuka (1968) and 
are shown in Chapter 2. A typical simulated profile is shown in Figure 5.6 
for a frequency f e of 5 KHz for loss tangent, the relative dielectric constant 






10,000 
















c '' f) simulation of typical lithospheric profile with thot in motel a ar-cvar chamber (From Iizuka.1969 ) 
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e re being 4 and constant with depth. The depth of the lithosphere D m was 
30 km, and a modelling scaling factor of 10^ was used so that the simulation 
chamber was 30 cm deep, D m = 0.3 meters. An aluminum plane represented the 
large loss-tangent upper crust, and a water solution of 0.2 mol NaCl repre¬ 
sented the mantle. The simulated profile of the conducting layer was 
obtained by diffusion of the NaCl, through a porous' film into the agar-agar 
at the bottom. 

Iizuka (1968, Chapter XX) describes how to fabricate profiles of 
various shapes including (a) complimentary error function, (b) error function, 
(c) linear profile, (d) periodic profile, and (e) miscellaneous profiles. 

5.3.2.1 Field strength measurements 

It should be mentioned that the experimental modelling chamber of Iizuka 
could be arranged also with non-parallel top and bottom, a vertical step 
profile (similar to Brown and Gangi, 1963), a vertical slant profile, or 
with simultaneous horizontal and vertical profiles. 

5.3.2.1.1 Non-parallel plate region with homogeneous medium 

The geometry is shown in Figure 5.7. Table I shows the dimensions and 
electrical properties of the medium for various geometers. Note that for 
loss tangents of 0.036, 0.23 and 0.88, the corresponding conductivities are 
0.073, 0.47 and 1.78 Si/m, the latter approaching that for sea water. 

The resulting electric field strength E in Figure 5.8 is shown plotted 
as E R vs the range R for a quarter-wave monopole radiator. Note maxima 


and minima for the lowest loss-tangent. The smaller the angle the smaller 






































5.19 


TA3L3 5.1 


Dimensions and electrical properties of the medium in a non-parallel plate 
region. 


At 114 MHz 


Coding 

h 

(in cm) 

Angle 
(in deg) 

* r 

tan 5 

3a 

8. 3 

3°53' 

78 

0. 036 

3b 

8. 3 

3°53' 

78 

0. 23 

3c 

8. 3 

3°53 ! 

78 

0. 88 

3a 

29. 0 

-3°53' 

78 

0. 036 

Tb 

29. 0 

-3°53' 

73 

0. 23 

3c 

29. 0 

-3°53' 

78 

0. 88 

5b 

7. 7 

-o 

0 

73 

0. 23 

5c 

7. 7 

5° 

73 

0. 83 

6a 

39. 0 

5°57' 

73 

0. 036 

6b 

39. 0 

5°57' 

78 

0. 23 

6c 

39. 0 

5°57 1 

79 

0. 88 
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FIG.5B PLOT OF E^R vs. R IN UNPARALLEL PLATE REGION 
FILLED WITH A CONDUCTING MEDIUM AT 114 MHz. 
SEE'THE PARAMETERS IN TABLE 5.1 (1.2 IWJ 
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the attenuation. As conductivity increases, the attenuation becomes less 
sensitive to angle changes. That is, at high conductivity the propagation 
depends more on the losses of the medium than by the geometry of the bottom. 
This is useful to know for operation at sloping shore lines. 

5.3.2.1.2 Field strength measurements in an inhomogeneous medium 

For the lithospheric guide model an agar-agar gel represented the 
lithosphere, covered with a metal plate top. A quarter-wave monopole 
extended through the plane into the medium. As such, the configuration 
was similar to a parallel waveguide with metallic walls except that the 
bottom was replaced by an inhomogeneously conducting wall. Figure 5.9 
is a representative result. The depth profile of e r and loss-tangent is 
shown in Figure 5.9 (b), while the amplitude and phase of the magnetic 
field strength are shown in Figure 5.9 (a). Note the model was 0.383 
X thick, yet the field strength decreases monotonically with depth. The 
field was measured at R = 50 cm; a perfect guide, less than half a wave¬ 
length thick would have a field magnitude invariant with depth. 

5.3.2.1.3 Comment - Images 

Figure 5.9 represents measurements on a model lithosphere for which 
the upper wall was a metal plate above the inhomogeneous lithospheric model. 
In essence, this could be considered one-half of a lithosphere having a 
symmetrical inhomogeneous profile, the upper half being the image of the 
lower one shown. Since the excitation was by a vertical monopo’e, the 
result would be that of a center-driven dipole, giving maximum field strength 


at the center where the conductivity (loss tangent) is minimum 
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Several of the more recently postulated profiles, such as those of 
Keller or of Housely discussed in Chapter 2, can be approximated by quasi- 
symmetrical profiles. Calculations of Field and Dore (1973) show maximum 
field strengths near that for minimum conductivity. 

Advantage could be taken of this image effect along with the various 
inhomogeneous model profiles achieved by Iizuka using the agar-agar gel 
with sodium chloride and ethyl alcohol diffusion discussed in Section 5.3.2. 
The uniform profiles of Brown and Gangi (1963) with saline water could 
likewise be treated - the scaling parameters would have to be reconsidered 
See Chapter 6. 

The use of images for measuring wall losses is also discussed in 
Chapter 6. Here, however, the question of similitude scaling application 
is difficult to answer because the walls are necessarily water instead of 
metal, in order to make distance or depth dependent measurements in the 
"walls" be the fields excited by dipoles in the lithosphere or in the 
walls. Nevertheless, the principles will be demonstrated and the simpler 
measurements understood before proceeding step by step to more complex 
cases. 

This leads to a discussion of modelling materials. 

5.4 Modelling Materials 

5.4.1 Typical High-conductivity Materials 

Fritschkneckt (1971) gives a table of typical high conductivity 
materials which might be considered for modelling purposes, and is shown 
in Table 5.2. Listed are values for metals, other solids, and some 
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Material 

Conductivity 
mho m 

Metals 


Aluminum 

1.34 ' 10 7 

Aluminum, commercial alloys 

1.22-3.60-10 

Brass 

1.2-1.3-1 O' 

Bronze, phosphor 

0.37-1.04> 10 7 

Bronze. “Everdur* 

J.3-5.3M0* 

Copper 

3.31 " 10’ 

Copper, tube (DHP copper* 

5.01' 10 7 

Lead 

4.53 ■ 10* 

Magnesium 

2.17 v 10 T 

M.ingantn 

2.27 ' 10’ 

Mercury 

1.04 ■ 10* 

''tool. nonmagnetic siumluss 

1.35-2.34 ' 10' 

Zinc 

1.02 • 10* 

Other solids 


Carbon 

2.0-2.3 • 10 4 

Graphite 

0.7- 1.2 v 10 5 

Kpow. commercial silver loaded 

5.0 v 10* 

Epovy. commercial copper loaded 

2.5 • 10 4 

-Lav-up* resin, graphite loaded 

0- 1.0 • '.0* 

’•Vast, bronze loaded 

-1.25- \0 4 

Silicone rubber, silver loaded 

1.3 • 10* 

L t qu ids 


\mmonium chloride. 


l'*% by weight in water it 13*C 

17.3 

<\»pper sulfate, is above 

J. J 

Hvdmohlonc acid, as above 

02.1) 

•vadium chlorine, ss above 

12.0 

Sulfuric acid, as above 

30.2 

Sea -.vator 

1-5 


5.2 Higher conductivity icdellins 
(fro*n Prischkneckt,19?1) 
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liquids. When large scaling ratios are required for models, a difficulty 
appears. One would like to have a material with conductivity between that 
for saline water and carbon, for certain frequency ranges. For distance- 
dependent field measurements, liquids or gels (such as agar-agar gels of 
Iizuka) are indicated. 

5.4.2 Low Dielectric Constant Solvents 

While searching for materials of low dielectric constants, Iizuka 
(1961) tabulated the results of his survey, shown here as Table 5.3. 

Dioxane was indicated but was too toxic and dangerous, so other solvents 
were used. The materials were desired in obtaining the performance of 
antennas in stratified media. The solvents were separated by thin films 
such as nylon (having low dielectric constants and conductivities). The 
horizontal stratification effects are reported by Iizuka (1961) and Iizuka 
and King (1962b). 

5.4.3 Higher Dielectric Constant Liquids 

For scaled modelling purposes,.higher dielectric constant liquids are 

desired. For this purpose, water is most often used, or ethyl alcohol or 

a mixture of the two, providing a range of 20 to 80 or so for relative 

dielectric constant. By adding salt, the conductivity can then be varied 
-4 

from 10 to about 10 Si/m. 

While less desirable for showing smooth profiles than the agar-agar 
mixtures of Iizuka, inhomogeneous profiles cam be simulated by horizontally 
stratified layers of solvents with differing electrical constants according 
to the profile desired. Such stratification might ultimately be considered 
in the schemes of Chapter 6 as one of the steps in building up from simpler 


models 
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6. Theory and Experiment In lithospheric Propagation: Models 
6.1 Introduction 

The possibility of wave-guide-like propagation of electromagnetic waves 
in the earth's lithosphere has been studied extensively by numerous investi¬ 
gators. The theoretical models are all quite similar and are based on ideali¬ 
zations of cross sections of the earth's crust such as that given by Levin 
[1966] and shown in Fig. 6.1. These idealizations in general include (a) a 
planar slab of a dielectric material with low conductivity to represent a 
lithosphere composed of granitic rock and basalt, (b) an upper layer of quite 
highly conducting material corresponding to the overburden on the continent 
and the ocean basin beyond this, and (c) a lower also highly conducting layer 
to represent the mantle. Refinements sometimes added include the air and the ion¬ 
osphere above the overburden and ocean as shown in Fig. 6.2, but their presence 
usually has a negligible effect. The excitation of electromagnetic waves in 
the dielectric slab and their reception at distant points has taken several 
forms. The first and most commonly described source consists of a vertical 
antenna in a borehole extending into the lithosphere or, short of this, ending 
in the overburden at some distance from the lithosphere. Wait's analysis [1966] 
shows that if the vertical antenna is not in the lithosphere or very close to 
it, adequate excitation is not possible. He suggests [1966, p.923] that "ever/ 
effort should be made to locate both the (vertical) source dipole and the obser¬ 
ver within the waveguide. Otherwise, some other mechanism of excitation should 
be used such as a horizontal electric dipole." This suggestion is followed by 
Frieman and Kroll [1973] whose Investigations include excitation not only by a 
vertical electric dipole but also by a horizontal electric dipole at the upper 
boundary of the lithosphere. They conclude that "the horizontal electric dipole 
case proves to be more interesting", that it "permits the utilization of much 
larger antennas and hence more efficient antenna design." They suggest that "a 
horizontal electric dipole is far superior to a vertical electric dipole." 








Fig.6.1 Profile of earth's crust according to S.B, Levin 


Ionosphere a 









The method of analysis of Frieman and Kroll is similar to that of Wait 
and others in that the relevant waveguide modes are determined. Their study of 
the propagating characteristics of the several modes led them to the conclusion 
that for a horizontal antenna and receiver on the ocean floor the behavior of 
the field "should be very similar to that of horizontal-dipole excitation near 
the plane boundary of two semi-infinite conducting media." However, it is well 
known [BaHos, 1966; King and Sandler, 1977] that transmission along such an 
interface is primarily by so-called lateral waves for which a complete theory is 
available. To date it has not been applied to an interface between two half¬ 
spaces with the properties of sea water and basalt or of the earth’s overburden 
and granitic rock. It appears, therefore, that for the idealized theoretical 
model consisting of a dielectric slab between conducting regions, analytical 
and numerical solutions of Maxwell's equations are available and should suffice 
to determine the nature of propagation between a pair of vertical or horizontal 
dipoles in the slab or near its boundary in the upper conducting layer. In 
particular, it should be possible to calculate the electric field at the re¬ 
ceiver as a function of distance, frequency, and the assumed electrical parameters 
of the media. Alternatives to such computations include,on the one hand, the 


very difficult and costly full-scale experiments and, on the other hand, suitably 
scaled and constructed experiments that satisfy all of the conditions of electro¬ 
magnetic similitude. Model measurements have the advantage over theoretical 
computations that they permit a variety of realistic modifications of the over¬ 
idealized theoretical models. 

In the absence of magnetic materials the conditions for modelling are: 


.22 ,2 2 

f e r ■ f e r 
m m m e e e 


far 2 
m m m 


far 2 
e e e 


( 6 . 1 ) 


where f is the frequency, r is a characteristic length, e - e 0 e r is the P erTIlit ' 
tivity and a is the conductivity. The subscript m denotes the model, the sub¬ 
script e the earth. Unfortunately such modelling is severely limited by the 


distances involved and by the materials and frequencies available for the 
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simulation and the experiments. A model of practical size requires a scale 

factor in the distances of the order of r /r ■ 10^ . With e - 4 to 8 for 

e m er 

the lithosphere and e ■ 80 for the water solutions most commonly used in the 

1/2 

model, 20 to 10. It follows that - (£ mr /£ er ) (r m /r e } " 

4 x 10“^ to 3 x 10~5 . Since accurate electrical measurements are difficult 

unless f = 1 GHz, it follows that f « (45 to 30) KHz. That is, completely 
m e 

scaled modeling of the lithosphere and its upper and lower bounding regions is 

limited to actual frequencies that satisfy f » 45 KHz. While this condition 

presents no practical limitation on VLF propagation in the lithosphere using 

vertical dipole excitation, it excludes all of the high frequencies which may 

be the most promising for lateral wave transmission using horisontal dipole 

excitation. The second condition for similitude requires that for e er ■ 4, a & /a m 

2 2 e. 

fr/fr »2x 10 If the lithosphere is modeled with water solutions as 
m m e e * 

-4 

was done by Brown and Gangi and by Iizuka, * 10 to 10 Si/m, so that 

o - 2 x 10" 1Q to 2 x 10“ 5 Si/m. This includes the entire useful range. On the 
e 

-2 -3 

other hand for an overburden with a * 10 to 10 Si/m, a m 3000 to 300 Si/m 

e m 

for which convenient materials are difficult to find. The same is true for 

the ocean with a - 4 Si/m and o = 4 x 10^ Si/m, with e__ * 80* 

e m er 

Actual models constructed by Brown and Gangi [1963] were designed to simu¬ 
late a frequency f * 2 to 4 KHz. The lithosphere was quite accurately scaled 

by salt water, the overburden or ocean and the mantle less accurately by metal 
or carbon plates. For waveguide propagation when the transmitting and receiving 
antennas are both vertical and embedded in the lithosphere, a less than exact 
scaling of the properties of the upper and lower conducting media may be adequate 
so long as the principal attenuation is due to losses in the lithosphere and not 
in the bounding walls. Note, however, that the use of solid plates for the over¬ 
burden or ocean is neither convenient nor sufficiently accurate when the exci¬ 
tation is by a vertical or horizontal dipole located at some point above the 
surface of the lithosphere. 
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The interesting models developed by Iizuka [1971,1972] have similar 
general limiations as those of Brown and Gangi with respect to the frequency 
range and the properties of materials. They are also restricted to excitation 
by dipoles in the model lithosphere. However, they do provide the means for 
producing a specified non-uniform conductivity profile across the lithosphere 
and for providing boundaries that are not plane parallel. In a restricted 
range of frequencies and with excitation by a vertical antenna within the 
lithosphere Iizuka's models offer a powerful tool for obtaining information 
about wave propagation in the earth insofar as the assumed profiles are good 
approximations of reality. Unfortunaely they are not readily generalized 
to permit the use of vertical or horizontal dipoles in the model overburden or 
ocean and for determining the electromagnetic field in this. 

6.2 Suggested Experimental Studies in Lithospheric Propagation 

Measurements on models combined with the principle of similitude provide 
a powerful and practical alternative to full-sized experiments in the earth 
and to theoretical studies. Owing to the very large scale factors involved 
and the limitations of experimental techniques and available materials the 
method is for all practical purposes limited to low frequency applications 
in the actual earth. Even for these the accurate modelling of the lithosphere 
virtually precludes a corresponding accuracy in the modelling of the overburden 
or the ocean and so requires the exciting and receiving antennas to be in the 
model lithosphere. In the actual earth this is possible only in boreholes 
through the overburden or into the ocean bed with vertical dipoles. In view of 
the theoretically derived conclusion by Frieman and Kroll [1973] that horizontal 
dipoles are a very promising alternative, future experimental studies should 
be concerned with a comparative investigation of the characteristics of both 
horizontal and vertical electric dipoles when located in the ocean near its 
interface with the lithosphere and not only with vertical dipoles in the litho- 
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sphere. Suggested experiments might include:a) vertical transmitting and 
receiving dipoles in the model lithosphere and at increasing heights in the 
model overburden or ocean; b) a vertical transmitting antenna in the model 
lithosphere and a horizontal receiving dipole in the model overburden or ocean 
near the interface with the lithosphere; and c) horizontal transmitting and 
and receiving antennas in the model overburden or ocean. In each case the 
receiving dipole is to be moved to measure the magnitude and phase of the 
electric field as the distance from the transmitting antenna is increased. 

Further measurements are indicated to determine the effect of non-parallel 
boundaries of the model lithosphere, boundary surfaces between the lithosphere 
and the overburden or ocean that are not smooth but undulate or are otherwise 
irregular or discontinuous, a layer of properly modelled sediment on the 
ocean bottom including the effect of embedding both transmitting and receiving 
antennas at various depths. 

Since communication between two horizontal dipoles in the overburden or 
ocean near its interface with the lithosphere is primarily by lateral waves 
[Banos, 1966; King and Sandler, 1977], experiments that involve such antennas 
must be designed to take maximum advantage of the somewhat unusual character¬ 
istics of this type of transmission. In particular, note must be taken of the 
fact that when the horizontal transmitting antenna is not far from the inter¬ 
face, the most useful component of the electric field is the radial one in a 
direction along the axis of the dipole. Furthermore, the optimum frequency 
for transmission is well above the low frequency range. Since accurate model¬ 
ling according to the principle of similitude cannot be carried out above full scale 
frequencies of two to three kilohertz, it is not practical to study lateral wave 
transmission in a comprehensive manner by models that satisfy the equations of 
similitude as formulated for the lithosphere. While this means that a completely 
scaled set of experiments is unavailable, it does not mean that meaningful and 
illuminating information cannot be obtained from suitably designed model experi- 


ments. These must, of course, be planned and interpretted in conjunction with 
the theory of lateral-wave transmission. 

In carrying out the suggested experiments it may be convenient and expedi¬ 
ent to separate the transmission losses in the model lithosphere from those in 
the model overburden, ocean, and mantle . This can be done initially by making use 
of a perfect dielectric (e.g. air or styrofoam) in the model lithosphere and sea 
water for the bounding regions. Subsequently a dielectric with losses and 
specified conductivity profiles comparable to those believed to obtain in the 
earth can be introduced in the simpler and more crucial experiments. 

6.3 Some Possible Laboratory Models for Lithospheric Experiments 

The simple model used by Brown and Gangi [1963] and the much more complicated 
one developed by Iizuka [1971,1972] stress the accurate modelling of the litho¬ 
sphere but represent the overburden or ocean and the mantle by highly conducting 
metal or carbon plates. In order to study the significance of the bounding 
regions in their effect on the transmission and simultaneously to make possible 
the embedding and moving of antennas in them, it is possible to model them with 
a liquid like water or alcohol or a mixture of these to obtain relative permit¬ 
tivities between 20 and 80 and to add salt to vary the conductivity in the 
-4 

range from 10 to 10 Si/m. The model lithosphere can be composed initially of 
air or styrofoam with negligible losses and subsequently by wood or alcohol with 
higher conductivity. The bounding styrofoam can be smooth, undulating, or 
arbitrarily irregular. Above the air or styrofoam is a metal plate that can be 
parallel to the surface of the water or inclined at an arbitrary angle. It 
serves as an image plane for the entire structure below it. This effectively 
provides a model lithosphere that has twice the thickness between the metal 
plate and the surface of the water. It is bounded above and below by salt 
water to model both the overburden or ocean and the mantle. 

A sketch of the suggested model is shown in Fig. 6.3. The vertical 
transmitting and receiving monopoles when combined with their images act as 
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Fig. 6.3 Suggested apparatus for lithospheric and lateral-wave 
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dipoles in the center of the model lithosphere. The receiving monopole travels 
in a slot to measure the amplitude and phase of the vertical electric field 
from near the transmitter to distant points. The upper plate can be inclined 
and the water-styrofoam surface can be made undulating ot otherwise irregular,Fig.6.4. 

Also shown in Fig. 6.3 are transmitting and receiving dipoles in the salt 
water. In the figure they are horizontal but in practice they can be oriented 
to be vertical or at any desired angle. The receiving element can be used to 
measure the amplitude and phase of the electric field in the model overburden 
or ocean. All components of the electric field can be determined. Since 
the dipoles have identical images which are symmetrically located in the model 


mantle, the field actually measured is due to both the source and the image. 
However, the contributions are readily separated theoretically, so that the 


field due to one source alone can be studied. 


In order to carry out the suggested measurements in a tank approximately 


3x2 meters, a frequency in the range from 0.3 to 1 GHz is required. This 
gives between three and ten wavelengths in air for studying the decay of the 
electric field in the model lithosphere when excited by a horizontal dipole 
in the overburden or ocean or by a vertical monopole in the model lithosphere. 
With the former the effect of varying the distance between the antenna and the 
bounding surface can also be investigated. The general arrangement shown in 
Fig. 6.3 provides the means for carrying out the suggested experiments and so 
obtaining direct information on the transmission of electromagnetic waves 
through the lithosphere and of lateral waves along its boundary under a vari¬ 
ety of circumstances arranged to simulate possible conditions in the earth. 
Such a program of experimentation and its correlation vzith theory would 
seem to be a logical and necessary continuation and extension of the basic 
work by Brown and Gangi and especially by Iizuka. 
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7. Summary and Conclusions 

Propagation in the earth's crustal waveguide has been studied by numerous 
investigators. The properties of available conceptual depth profiles of elec¬ 
tric constants, mainly conductivity, are reviewed. The propagating region was 
assumed to be a planar slab, bounded on top by an ocean or overburden and on 
the bottom by the mantle or Moho, the bounding regions being of very large 
loss-tangent. 

The waveguide theories extant are then summarized. Earlier analyses 
assumed a slab of uniform depth profile. Many of the investigators illustrate 
the modal properties for modes of low attenuation, others give samples of 
field calculations. Usually the exciting sources were vertical electric 
dipoles (VED) extending through bore-holes in the overburden (or ocean) into 
the lower loss-tangent crustal waveguide. Analyses for non-homogeneous pro¬ 
files included the modal solutions of Wait for an exponential conductivity 
profile, and those of Schwering, Peterson, and Levin for an inverse-square 
depth profile of both dielectric constant and conductivity. Full wave modal 
solutions were employed by Field and Dore using recent conceptual ”bell-shaped" 
profiles of Keller and of Housley which showed minima and very much smaller 
conductivities at depths of the order of 10 km. 

Peking VED excitation in the waveguide, Wait suggested HED excitation, 
a suggestion which was followed up by Frieman and Kroli who showed that for 
ocean emplaced antennas at the ocean-lithosphere interface, the HED was, 
superior to the VED. Excitation should be like that for HED excitation near 
the plane boundary of two semi-infinite media. For dipoles in the ocean 
near the ocean-air boundary, King and Sandler extended the work of Banos. 
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Propagation is by means of lateral waves. This work was extended to the case 
where the air was replaced by the lithosphere. Calculations were made for a 
wide range of parameters for the radial component of electric field from the 
HED, which component is the useful one especially for intermediate and large 
ranges. The mechanism may be described as "down-under-and-up" in consonance 
with the popular notation of "up-over-and-down" for HED sources in the sea 
near the sea-air boundary. The advantages and disadvantages of the lateral 
wave transmission were discussed for the ocean-lithosphere case. 

The relations for the scaling of models according to electrodynamic 
similitude are presented. Examples are given from the work of Brown and Gangi 
for homogeneous depth profiles where the modelled crustal slab (saline) water 
is bounded on top by metal plates. For non-homogeneous profiles, examples 
are given from the work of Iizuka who used an aaar-agar gel with diffusion of 
sodium chloride to vary the loss-tangent with depth. The top boundary was a 
metal plate, the bottom being the higher conductivity sodium chloride diffuser 
Measurements were made with planar and sloping boundaries. 

To study the effects of bounding regions on transmission including losses 
and simultaneously make measurements of the effects of embedding and moving 
antennas in them, the laboratory models should have "walls" composed of a 
liquid like water or alcohol or a mixture of them which can vary the relative 
dielectric constant from 20 to 80 and to add salt to vary the conductivity 
from 10 to 10 Si/m. Above such a wall would be the lithosphere, initially 
air or styrofoam, and subsequently a substance like wood or alcohol to in¬ 
crease the complex dielectric constant. Above the "lithosphere" is a metal 
plate which can be parallel or inclined. The styrofoam boundary can be 
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undulating or irregular. The metal plate, which has movable probes, acts 
as an image plate. 

It is recommended that an orderly build-up of modelling experiments be 
initiated with a laboratory model just described, and continuing with more 
complicated depth profiles. Although available materials may limit accurate 
modelling according to the principles of similitude, meaningful and illumi¬ 
nating information can be obtained for a better understanding of lithospheric 
transmission including lateral waves and losses. 

The more complicated crustal depth profiles require a much more elaborate 
set-up possibly patterned after that of Iituka. For measurements and their 
interpretation then, there will be a need for scientists of considerable 
experience including support personnel. The program should be long range 
in anticipation of the development of more accurate depth profiles from full 
scale geophysical experiments and analyses. 
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